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tj.  AfSfRACT  (MttmMim 200woirMi  The  Biot  theory  of  propagation  in  a  porous  medium  provides  a 
mathematical  framwork  for  studying  acoustic  interaction  with  the  seafloor.  The 
theory  considers  the  two-phase  porous  nature  of  marine  sediments  in  contrast  to  the 
classical  models  of  wave  propagation  in  the  seafloor  that  consider  marine  sediments 
as  an  extended  single-phase  fluid  or  solid.  A  boundary  value  problem  is  set  up  and 
solved  for  a  line  source  in  a  fluid  medium  above  a  poro-vicroelastic  halfspace. 
Expressions  for  the  reflected  ^d  transmitted  field  are  given  in  integral  form  and 
asymptotic  expansions  in  the  high-frequency,  far-field  limit.  A  set  of  simultaneous 
equations  is  solved  to  give  plane  wave  reflection  and  transmission  (Type  I,  Type  II 
and  shear  wave)  coefficients.  These  equations  also  yield  the  Scholte,  pseudo-Scholte 
and  pseudo-Rayleigh  wave  phase  velocities  and  attenuations.  The  plan  wave  coefficieni  s 
and  the  surface  wave  velocities  and  attenuations  are  compared  with  commensurate 
quantities  from  the  single-phase  theories.  A  number  of  recent  experiments  of  high 
frequency  transmission  through  a  water-sand  interface  have  indicated  anomalously  high 
transmitted  energy  at  angles  near  the  critical  angle.  The  Type  IT  wave  (predicted  by 
Biot  theory  but  not  the  single-phase  theories)  was  suspected  as  a  possible  reason  for 
the  anomalies.  Data  from  one  of  the  experiments  are  examined  using  Biot  theory  to 


(cont) 


determine  whether  or  not  the  Type  II  wave  was  responsible  for  the  anomalies. 
Calculations  indicate  that  the  contributions  from  the  Type  II  wave  are 
negligible  and  thus  not  the  cause  of  anomalies.  Much  of  the  anomalous  behavior 
appeared  to  arise  simply  as  a  result  of  an  incorrect  compressional  wave  attenuation 
in  the  geoacoustic  model. 


ABSTRACT 


ill 


The  Biot  theory  of  propagation  in  a  porous  medium  provides  a  mathematical 
framework  for  studying  acoustic  interaction  with  the  seafloor.  The  theory  considers  the 
two-phase  porous  nature  of  marine  sediments  in  contrast  to  the  classical  models  of  wave 
propagation  in  the  seafloor  that  consider  marine  sediments  as  an  extended  single-phase 
fluid  or  solid.  A  boundary  value  problem  is  set  up  and  solved  for  a  line  source  in  a  fluid 
medium  above  a  poro-viscoelastic  halfspace.  Expressions  for  the  reflected  and  transmitted 
field  are  given  in  integral  form  and  asymptotic  expansions  in  the  high-frequency,  far-field 
limit.  A  set  of  simultaneous  equations  is  solved  to  give  plane  wave  reflection  and 
transmission  (Type  I ,  Type  n  and  shear  wave)  coefficients.  These  equations  also  yield  the 
Scholte,  pseudo-Scholte,  and  pseudo-Rayleigh  wave  phase  velocities  and  attenuations. 
The  plane  wave  coefficients  and  the  surface  wave  velocities  and  attenuations  are  compared 
with  commensurate  quantities  fiom  the  single-phase  theories. 

A  number  of  recent  experiments  of  high  frequency  transmission  through  a  water- 
sand  interface  have  indicated  anomalously  high  transmitted  energy  at  angles  near  the  critical 
angle.  The  Type  n  wave  (predicted  by  Biot  theory  but  not  the  single-phase  theories)  was 
suspected  as  a  possible  reason  for  the  anomalies.  Data  from  one  of  the  experiments  are 
examined  using  Biot  theory  to  determine  whether  or  not  the  Type  n  wave  was  responsible 
for  the  anomalies.  Calculations  indicate  that  the  contributions  from  the  Type  n  wave  are 
negligible  and  thus  not  the  cause  of  the  anomalies.  Much  of  the  anomalous  behavior 
appeared  to  arise  simply  as  a  result  of  an  incorrect  compressional  wave  attenuation  in  the 
geoacoustic  model. 
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Chapter  1 
INTRODUCTION 

An  important  problem  in  the  field  of  underwater  acoustics  is  the  interaction  of  a 
sound  field  with  the  seafloor.  The  seafloor  is  an  acoustically  complex  medium  exhibiting 
properties  that  can  be  highly  variable  in  both  vertical  and  lateral  extent  On  a  microscopic 
scale  the  seafloor  is  no  less  complicated;  marine  sediments  manifest  themselves  as  two-  and 
sometimes  three-phase  materials  with  solid  grains  of  variable  dimension,  shape,  and 
composition  randomly  oriented  and  connected  in  space,  the  interstices  being  filled  by 
seawater  and  dissolved  gases. 

1.1  Background 

The  modeling  of  wave  propagation  through  the  physically  complex  medium  of  the 
seafloor  has  progressed  from  treating  the  sediments  as  a  lossy  fluid,  to  a  viscoelastic  solid, 
to  a  poro- viscoelastic  solid.  The  great  bulk  of  modeling  efforts  of  acoustic  interaction  with 
the  seafloor  in  the  ocean  acoustics  community  have  considered  the  seafloor  sediments  to  be 
modelable  as  a  fluid  [e.g.,  Ingenito  (1973),  Ellis  and  Chapman  (1980),  Rubano  (1980), 
Dicus  and  Anderson  (1982),  and  Brocher  (1983)].  The  fluid  model  has  been  sufficient  in 
fact  to  successfully  predict  various  measures  of  a  propagating  acoustic  field:  most 
commonly  transmission  loss.  Some  of  the  experimental  results  [e.g.,  Ingenito  and  Wolf 
(1976),  McDaniel  and  Beebe  (1980),  Chapman  and  Ellis  (1984),  Beebe  and  Holland 
(1986),  Hughes  et  al.  (1990)],  however,  indicated  that  the  effect  of  the  sediment  rigidity 
could  be  significant  for  various  environments.  Modeling  efforts  [Fryer  (1978),  Vidmar 
(1980),  Holland  (1985)]  that  considered  the  sediments  as  viscoelastic  attempted  to  determine 
the  bounds  of  validity  for  the  fluid  approximation  as  a  function  of  frequency  and  seafloor 
environment 
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In  the  meantime  a  substantial  effort  was  directed  toward  an  even  more  rigorous 
accounting  of  the  physics  of  propagation  in  marine  sediments,  using  Biot  theory  (1956a, 
1956b,  1962),  which  considered  marine  sediments  as  a  poro-elastic  solid.  The  theory 
predicts,  in  contrast  to  both  the  fluid  model  and  the  viscoelastic  model  (which  will  be 
referred  to  as  single-phase  models)  a  non-linear  frequency  dependence  of  wave  attenuations 
and  a  compressional  wave  of  the  second  kind  (which  will  hereafter  be  referred  to  as  the 
Type  n  wave)  in  addition  to  the  classical  compressional  (Type  I)  and  shear  wave  of 
geophysics. 

The  foundations  of  the  Biot  theory  were  laid  by  Rayleigh  [Strutt  (1883)]  who 
considered  propagation  in  a  porous  medium  whose  solid  frame  was  rigid  and  whose 
Interstices  were  circular  cylinders.  Later,  Zwikker  and  Kosten  (1949)  developed  the  rigid 
frame  theory  further  and  began  initial  developments  of  a  flexible  frame  theory.  It  was  Biot, 
however,  who  provided  the  first  comprehensive  theory  of  wave  propagation  in  fluid 
saturated  porous  media  that  accounted  frar  both  the  losses  due  to  the  relative  motion  between 
the  interstitial  pore  fluid  and  the  solid  frame  and  accounted  for  displacements  in  the  solid 
frame. 

In  an  exhaustive  series  of  papers,  Deresewicz  and  co-authors  (1962-1964)  used 
Biot  theory  to  examine  a  number  of  seafloor  interaction  problems.  Later,  Stoll  (1974) 
introduced  the  notion  of  complex  frame  moduli  to  account  for  viscoelastic  effects  and  in  a 
series  of  papers  (Stoll  (1977),  Stoll  (1979),  Stoll  (1980),  Stoll  and  Kan  (1981)]  was 
largely  responsible  for  bringing  Biot  theory  to  the  attention  of  the  underwater  acoustics 
community.  Thus,  the  poro-viscoelastic  model  (sometimes  termed  the  Biot-Stoll  model) 
considers  loss  mechanisms  due  to  the  relative  motion  between  the  fluid  and  frame  as  well 
as  (approximately)  constant  Q  effects  such  as  the  squeeze-film  phenomena,  frictional  losses 
at  grain-to-grain  contacts,  and  relaxation  phenomena  associated  with  electrochemical 
effects.  Yamamoto  (1983a,  1983b)  was  also  active  in  applying  the  theory  to  a  number  of 
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problems  and  generalized  the  form  of  the  visco-dynamic  operator  to  include  a  distribution 
of  grain  sizes  [Yamamoto  and  Turgut  (1988)]. 

Laboratory  experiments  have  to  date  confirmed  many  aspects  of  the  Biot  theory 
including  the  fi^equency  depenf^-nce  of  the  Type  I  comprcssional  velocity  and  attenuation 
[e.g.,  Hovem  and  Ingram  (1979),  Addison  (1984),  Winkler  (1985)],  the  dependence  of  the 
Type  I  compressional  wave  velocity  and  attenuation  on  the  properties  of  the  interstitial  pore 
fluid  [e.g.,  Domenico  (1977)  and  Costlcy  (1986)],  the  frequency  dependence  of  the  shear 
wave  attenuation  [Brunson  and  Johnson  (1980)]  and  the  existence  of  the  Type  n  wave 
[Paterson  (1956)  and  Fiona  (1980)].  Actually,  the  existence  of  the  Type  n  wave  was  well 
known  in  aeroacoustics  and  its  properties  had  been  measured  in  porous  architectural 
materials  such  as  foams  and  fibrous  tiles.  At-sea  experiments  [e.g.,  Ingenito  (1973)]  also 
have  indicated  that  the  firequency  dependence  of  the  Type  I  wave  compressional  wave 
attenuation  is  not  always  linear  with  frequency  and  Biot  theory  has  been  successfully 
employed  in  explaining  such  non-linearities  [Beebe  (1981)].  An  excellent  recent  summary 
of  the  state  of  understanding  of  the  frequency  dependence  of  compressional  wave 
attenuation  in  marine  sediments  and  the  evidence  for  a  non-linear  fiiequency  dependence  is 
given  by  Kibblewhite  (1989). 

1.2  Objective  and  Approach 

The  motivation  for  this  work  was  to  understand  how  the  porous,  two-phase  nature 
of  marine  sediments  influences  reflection  from  and  transmission  through  the  seafloor  since 
the  majority  of  current  models  of  seafloor  interaction  consider  marine  sediments  to  be  a 
single-phase  fluid  or  a  single-phase  solid.  In  order  to  focus  on  the  effects  of  the  two-phase 
poro-viscoelastic  medium,  a  great  many  complexities  in  nature  are  ignored  here,  including 
interface  roughness  and  slopes;  velocity,  attenuation,  and  density  gradients;  and 
macroscopic  inhomogeneities. 
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The  approach  taken  is  to  set  up  and  solve  a  boundary  value  problem  of  a  fluid 
(seawater)  overlying  a  poro-viscoelastic  halfspace  (sediment)  where  the  source  is  in  the 
upper  medium.  For  simplicity  the  source  is  chosen  to  be  a  line  source  such  that  the 
problem  can  be  solved  in  two  dimensions.  The  linearized  wave  equation  is  employed  in 
the  upper  halfspace  and  the  equations  of  motion  for  a  poro-viscoelastic  medium  derived  by 
Biot  are  employed  for  the  sediment  The  form  of  the  Biot  theory  in  this  work  largely  follows 
Stoll  (1974).  The  effects  arising  from  the  two-phase  natme  of  the  sediments  in  the 
reflection/transmission  problem  are  isolated  by  a  comparison  with  commensurate  results 
from  the  single-phase  models.  The  comparison  of  reflection  and  transmission  between  the 
single  phase  models  and  the  poro-viscoelastic  model  is  performed  for  three  sediment  types: 
a  silty  clay,  a  carbonate  sand,  and  a  sandstone.  These  sediment  types  were  chosen  as 
representative  of  common  sediment  types  and  because  they  represent  distinct  relationships 
between  the  ratio  of  their  velocities  to  that  of  the  overlying  bottom  water.  For  the  silty  clay, 
the  Type  I  and  shear  wave  velocity  are  both  less  than  that  of  the  upper  halfspace;  for  the 
carbonate  sand,  the  Type  I  wave  velocity  is  greater  and  the  shear  velocity  less  than  that  of 
the  upper  halfspace;  and  for  the  sandstone,  the  Type  I  wave  and  shear  velocity  are  both 
greater  than  that  of  the  upper  halfspace.  The  Type  n  velocity  is  less  than  that  of  the  upper 
halfspace  for  all  sediment  types. 

In  Chapter  2,  the  time  harmonic  equations  of  motion  are  written  for  both  media  and 
Fourier  transformed  into  wavenumber  space.  The  dispersion  relations  for  the  poro- 
viscoelastic  medium  are  obtained  and  examples  of  the  frequency  dependent  velocities  and 
attenuations  for  the  three  sediment  types  are  shown.  The  boundary  conditions  are  then 
applied,  resulting  in  a  set  of  simultaneous  equations  that  can  be  solved  for  the  plane  wave 
reflection  and  transmission  coefficients.  In  Chapter  3,  the  reflection  coefficients  from  a 
poro-viscoelastic  halfspace,  a  viscoelastic  halfspace,  and  a  fluid  halfspace  are  compared  in 
frequency  and  angle  space  for  the  three  kinds  of  sediments.  Differences  in  the  reflectivity 
between  models  is  noted  and  observations  are  made  regarding  tire  effect  of  using  one  model 
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or  another  in  acoustic  modeling.  Several  approximate  models  are  also  examined.  The 
secular  equation  is  derived  and  surface  wave  speeds  and  attenuations  are  obtained  and 
compared  with  predictions  from  the  single-phase  viscoelastic  solid  sediment  model.  In 
Chapter  4,  the  inverse  transform  is  applied  using  a  combination  of  techniques  to  obtain  the 
contribution  from  each  wave  type  (i.e.,  direct  path,  geometrical  reflected/transmitted  ray, 
head  wave,  surface  wave)  and  the  results  are  compared  with  the  “exact”  solution  (obtained 
by  numerical  integration)  and  with  the  solution  obtained  by  a  Fast  Fourier  Transform.  The 
decomposition  of  the  field  into  discrete  wave  types  allows  identific  ation  of  dominant  path 
contributions  that  arise  in  specific  source/receiver  geometries.  In  Chapter  5,  recent 
anomalous  experimental  data  of  acoustic  transmission  through  a  water-sand  interface  are 
examined  using  Biot  theory  to  determine  whether  the  anomalies  are  attributable  to  effects 
predicted  by  Biot  theory.  A  summary  and  conclusions  of  this  work  are  given  in  Chapter  6. 


Chapter  2 

THEORY  I  -  FORWARD  TRANSFORM 
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In  this  chapter  the  equations  of  motion  are  given  for  each  medium  and  transformed 
to  the  frequency  domain,  the  boundary  conditions  are  stated  and  a  set  of  simultaneous 
equations  are  written  to  obtain  the  reflection  and  transmission  coefficients. 

2.1  The  Fluid 

The  linearized  equation  of  motion  for  an  ideal  fluid  with  a  line  source  perpendicular 
to  the  x-z  plane  can  be  written: 

V2p  -  ^  ^  =  Po  8(x)  8(z-h)  1^  f(0  (2.1) 

where: 

p  is  the  dynamic  acoustic  fluctuation  pressure 
Vg  is  the  acoustic  wave  speed 

Po  is  the  static  density 

f(t)  is  the  forcing  function 

and  the  (line)  source  is  located  at  x=0,  z=h.  Letting  f(t)  be  a  unit  time  harmonic  function, 
f(t)  =  e  equation  (2.1)  becomes 

V2p  +  kg^p  =  -©2p^  5(x)  5(z-h)  ;  ko  =  (2.2) 

’O 

where  the  harmonic  time  dependence  has  been  suppressed.  Defining  a  Fourier  Transform 
on  the  X  variable: 

- -  oo 

p(kx,z)=  Jp(x,z)e*»*dx 
— 

p(x,z)  =  ^  Jp(k^.z)e~‘^^^dk^ 


(2.3) 


then 
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(2.4) 

(2.5) 

(2.6) 

(2.7) 


The  integrand  has  simple  poles  at 

•^  =  ±(ko2-kx2)l/2  (2.8) 

Taking  the  case  (z-h)  <  0  first  and  closing  the  contour  in  the  upper  half  plane 


Figure  2. 1  Path  of  integration  in  the  complex  kz  plane  for  (z-h)  <  0 


8 


we  have  from  the  Residue  Theorem  that 


(2.9) 


(2.10) 


To  make  the  integral  on  Cr  vanish  as  R-4  «»;  we  require  that  Im  (kz)  >  0  or, 
Im{(k:o^-  kx^)''^^<  0 ;  then 


p(kx,z<h)  = 


2 


-i(k2-k2)l/2(2_h) 


(2.11) 


For  (z-h)  >  0:  we  close  the  contour  in  the  lower  half  plane 


Figure  2.2  Path  of  integration  in  the  complex  kz  plane  for  (z-h)  >  0 
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and  from  the  Residue  Theorem 


R 

J 

-R 


^-.k^(z-h) 

k^  +  k^-k^ 

Z  X  o 


dk  + 

z 


^-ik^(z-h) 

k^  +  k^-k^ 

Z  X  o 


dk  =  — 2;r£res 

Z 


hJ 


(2-12) 


(2.13) 


Note  that  the  minus  sign  on  the  RHS  of  Eq.  (2.12)  follows  from  Cauchy’s  law  closing  the 
contour  clockwise.  To  make  the  integral  on  Cr  vanish  as  R->0  we  require  that  Im{kz}<  0 
orim  {(ko^  -  kx^  then 


^  ux  10)  Po  c 

p(kx,z^h)  = - ^ 

2  (kj-k^) 


W 


(2.14) 


and  for  all  z  (combining  Eq.  (2. 1 1)  and  Eq.  (2.14)) 


.  2  -i(k?-lcj)*^2|z-h| 

,  »0)^PoC  '  o  ' 

p(k„z) - 


(2.15) 
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The  total  solution  is  the  sum  of  the  homogeneous  and  the  particular  solution.  The 
radiation  condition,  however,  requires  that  Im  {(ko^  -  <  0  so  that  Bi  =  0,  and  the 

full  solution  is 


p(k^,z)  =  Aie  ^0 


Now  let 


Ai=9l 


tpoCO^ 


.2(ko2  -  kx2)l/2J 


(2.16) 


where  SI  can  be  thought  of  as  a  reflection  factor  for  the  amplitude  of  the  velocity  potential, 
then 


•*/i_  \  Po 


+  g-‘(k2-k2)>/2|z-h| 


(2.17) 


The  boundary  conditions  will  require  an  expression  for  the  field  in  terms  of  th(! 
displacement  vector  Uq  where  the  subscript  denotes  the  upper  medium.  From  Eulers 
equation 

-Vp  =  Po^  (2.18) 

Since  we  arc  assuming  a  time  harmonic  field  (e  »“*) 


Vp  =  poto2Uo 


(2.19) 
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and  transforming  Euler’s  equation,  we  get 


«  dp  2 
-a,tkxP+a3  — =  po0)  Uo 


where  a i,  33  are  the  x  and  z  unit  vectors  respectively.  Letting  a  =  (ko^  - 

f  ..X 


(2.20) 


2a 


a3fli^j[9i(-ia)e  -'“"±  (-ia)c  } 


(2.21) 


Note  that 


iLg-ialz-hl  =.j-ote-»alz-hl  f^iz-bl 
dz 


and 


+ 1  for  z  >  h 
- 1  for  z  <  h 


For  the  last  term  of  Eq.  (2.21)  the  minus  sign  is  chosen  because  we  want  to  use  this 
expression  for  the  boundary  conditions  (where  z  =  0)  so  that 


(2.22) 
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2.2  The  Poro-viscoelastic  Halfspace 
2.2.1  Equation  of  Motion 

In  this  section  the  equations  of  motion  for  a  pxjro-viscoelastic  medium  are  derived, 
the  dispersion  relation  is  found  and  tne  Fourier  transformed  displacement  vectors  for  the 
motion  of  the  sediment  frame  and  the  relative  motion  between  the  frame  and  the  interstitial 
fluid  are  obtained. 

The  equations  of  motion  are  obtained  by  specifying  both  the  constitutive  equations 
(relating  the  internal  restoring  forces  of  the  medium  to  the  corresponding  deformations)  and 
the  force  balance  equadons  which  come  from  Newton’s  second  law.  For  the  two-phase 
material  under  consideration  each  kind  of  conservation  law  is  expressed  by  two  equations; 
one  for  the  solid  component  and  one  for  the  fluid  component. 

For  the  constitutive  equations,  generalized  Hooke’s  law  for  an  isotropic  porous 
medium  yields  in  dyadic  notation 

X  =  I  {>.(V  u)  -  C  (V  w))  (V  u  -m  V)  (2.23) 

where  T  is  the  stress,  I  is  the  identity  matrix,  X  and  |i  are  the  Lame’  constants,  u  is  the 
absolute  displacement  of  the  solid  frame,  w  is  the  relative  displacement  of  the  two  phases 

w  =  p  (u  -  U)  (2.24) 

U  is  the  absolute  displacement  of  *he  fluid,  P  is  the  material  porosity,  and  C  is  an  operator 
that  characterizes  the  elastic  and  inelastic  response  of  the  sediment  frame  (see  Appendix  A). 
Without  the  term  C(V  w),  Eq.  (2.23)  reduces  to  the  constitutive  equation  for  an  isotropic 
elastic  solid.  The  constitutive  equation  for  the  fluid  is: 


Pf  =  M(V  w)  -  C(V  u) 


(2.25) 


where  pf  is  the  pressure  in  the  interstitial  fluid  and  M  is  an  operator  of  the  same  form 
described  for  C.  Eq.  (2.25  )  has  the  same  form  as  the  equation  of  continuity  for  an 
isotropic  fluid. 


The  force  balance  equations  likewise  resemble  their  isotropic  counterparts, 

Vx  = 

-Vp  = 

where 

pT  =  the  total  static  density  of  the  medium 

Pf  =  the  density  of  the  interstitial  fluid 

Pb  =  the  complex  density 

The  total  static  density  of  the  medium  can  be  obtained  by  a  simple  material  average,  i.e. 

PT  =  pfp  +  ps(l-  p)  (2.28) 

where  ps  is  the  density  of  the  solid  grains.  The  complex  density  was  derived  by  Biot 
(1962)  as 

Pb  =  c^-ia)^F(Y)  (2.29) 

The  first  term  accounts  for  the  fact  that  since  the  pore  space  is  tortuous  not  all  the  fluid 
moves  in  the  direction  of  the  applied  pressure  gradient.  The  parameter  c  is  equal  to  one 
where  the  pores  are  aligned  in  the  direction  of  the  pressure  gradient  and  greater  than  one  for 
randomly  ordered  pore  geometries  characteristic  of  sediment  interstices.  The  second  term 
accounts  for  viscous  drag  caused  by  fluid  viscosity  q  scaled  by  the  permeability  The 


~  [PtU  -  pfw]  (2.26) 

H  [pfu  -  Pbw]  (2.27) 
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function  F(y)  was  derived  originally  by  Zwikker  and  Kosten  (1949)  to  account  for  the 
transition  between  Poiseuille  flow  at  low  frequencies  and  Helmholtz  flow  at  high 
frequencies.  The  function  is  defined  as: 


where 


YT(y) 

1-2T(y)/Y 


(2.30) 


_  ber^(Y)  +  i  bei^(Y) 
^  ber(Y)  +  i  bei'(Y) 


(2.31) 


and 

Y=a((opf/Ti)l/2  (2.32) 


and  a  is  a  parameter  with  the  dimension  of  length  that  depends  on  the  size  and  shape  of  the 
pores. 

Now,  substituting  the  constitutive  equations  into  the  force  balance  equations 

V  •  t  =  V  {  A,  (V  u)  -  C(V- w) }  +  p.V2  u  +  |iV  (uV)  (2.33) 

and  since 

V  (uV)  =  V  (V  u)  (2.34) 

V  •  X  =  (H-  p)  V  (V  u)  -  CV(V  w)  +  pV^u  =  -o)2|'p.j.u  -pfw]  (2.35) 

and 

V  p  =  CV(V  u)  -  MV(V  w)  =-  ci)2[pfu  -Pbw]  (2.36) 

where  H  =  A.  +  2  p  . 

We  wish  to  express  the  equations  in  terms  of  potentials,  and  performing  the  usual 


Helmholtz  decomposition,  let 
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u  =  V<l)s  +Vx\|rs 


V-Vs  =0 


(2.37) 


then 


w  =  V4)f  +  Vx\jff  ;  =  0 

(H-n)  V  (V2<()s )  -  CV(V2<j)f)  +  V2(V<t»s+  Vxvs)  = 


(2.38) 


-0)2  [pt(V<1>s+  VxVs)  -  pf(V4)f  +  Vxvf )] 


CV(V2(j)s)  -  MV(V24)f)  =  -0)2  [pf  (V(l)s+  VxYs)  -  pf(V(j>f+  Vx\|rf)]  (2.39) 

The  irrotational  field  can  be  obtained  by  taking  the  divergence  of  Eqs.  (2.38, 2.39), 
i.e. 

HV2<|>s  -  CV24>f  =  .©2  Pf  <j,f]  (2.40) 

CV2<j)s  -  MV2{j)f  =  -0)2  [pf  <j)5-  pb<^f]  (2.41) 

and  the  rotational  component  of  the  field  can  be  obtained  by  taking  the  curl  of  Eqs.  (2.38, 
2.39) 

^L  V2i|fs  =-  0)2  [ptVs  -  PfVf]  (2.42) 

PfVs  =  PbVf  =  [c^  -  ^  Vf  (2.43) 

The  equations  for  the  irrotational  component  of  the  field  allow  solutions  of  the  form 

<l)s(r)  =  Ae-KqT)  (2.44) 

<}>f(r)  =  Be'KqT) 


where  q  •  q  =  k2 
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and  plugging  them  into  Eqs.  (2.40, 2.41)  we  obtain  the  dispersion  relation 

-k^H  +  CD^  Pt  -  0)2  pf 

-k^C  +  0)2  pf  k2M  -  0)2  pb 

or  in  terms  of  the  velocity  V=o)/k  we  have 


=  0 


(2.45) 


(Pf^  *  PTPb)  +  [PbH  +  PtM  -  2pfC]  +  C!2  -  HM  =  0  (2.46) 

There  are  two  solutions  to  Eq.(  2.46)  corresponding  to  two  dilatational  waves  (Type  I  and 
Type  II)  that  can  propagate  in  porous  media.  The  Type  I  wave  is  the  wave  of  classical 
geophysics  with  wave  speeds  on  the  order  of  .95  -  1.3  times  the  speed  of  the  interstitial 
fluid  speed  for  unconsolidated  marine  sediments.  Henceforth  the  subscript  to  designate 
this  wave  type  will  be  "p"  as  a  reminder  that  it  is  the  classical  p-wave.  The  Type  n  wave  is 
diffusive  (that  is  the  imaginary  part  of  the  wavenumber  is  not  small  compared  with  the  real 
part)  at  low  frequencies  for  water  saturated  sediments  and  propagatory  at  high  frequencies. 
For  this  wave  type  the  interstitial  fluid  and  the  sediment  move  essentially  out  of  phase.  The 
subscript  designating  this  wave  type  will  be  "2".  The  subscript  for  the  shear  wave  will  be 
"s".  Thus,  the  velocities  for  the  Type  I,  Type  II  and  shear  wave  will  be  Vp,  V2  and  Vs 
respectively.  Rewriting  the  potentials 

<t)s  (r)  =  Ape'^^^pT)  +  A2e’‘  (92‘'‘)  (2.47) 

<l>f(r)  =  Bpe-<(qpT)+  B2e-i(q2T) 

we  can  solve  for  the  ratio  of  coefficients  between  the  two  potentials  using  Eq.  (2.45). 
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Gps2a  =  PT 

Ap  C/Vp2  -  Pf 


(2.48) 


-§2  _  H/V2^  -pT 

^  A2  C/V22  -  pf 


(2.49) 


The  dispersion  relation  for  the  rotational  component  of  the  field  may  be  had  by  assuming 
solutions  of  the  form 


Vs  =  De'^Q*") 


(2.50) 


Vf  =  Ee'^Q’*") 


so  that  Eqs.  (  2.42,  2.43)  yield 


-k2p.  +  <b2  p.j.  .  0j2  Pf 


(2.51) 


The  shear  velocity  Vs  =  oVk  is  : 


Vs  = 


—  fv  2/ 


Pr-Rf 


(2.52) 


Gs  =  S  = 


(2.53) 


We  can  write  the  transformed  displacement  potentials  now,  letting  the  qi  have  z 
components  only 


qi  =  -33  (ki^  -  kx2)  1/2  =  -33  ai  t  =  p,  2,  s 


(2.54) 
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for  i  =  Type  1, 11,  and  shear  waves.  The  x  component  of  the  wave  vector  is  the  same  in 
both  the  fluid  and  porous  solid  (Snell’s  law)  and  we  choose  aito  be  negative  for  waves 

traveling  in  the  negative  z  direction  (away  &om  the  boundary),  so  that 
(j)s  =  Ape  ‘  “p^  +  A2e‘ 

<j)f=ApGpe‘“pz  +  A2G2e‘«2z  (2.55) 

Vs  =  a2De‘  “s^ 


Vf=  a2DGse‘®s^ 
Recalling  Eq.  (2.37),  then 


u(x,y)  =  ai  ^  +  as  ^  +  83 


ai  ^ 

X  *  dz 


w(x,y)=ai  +  83  ^  +  a3  ai 


and  from  our  definition  of  the  Fourier  Transform  [Eq.  (2.3)] 


I  ^ 


-  I  1  1  3  Vs)  I 

u=ail-«kA— ^1+83  -*xVs+-^ 


w  =  ai 


3Vf 

dz 


+  83 


-^xVf+ 


dz 


so  that 


u=-  t[ai{k;^Ape‘“p^ +kxA2e‘“2^  +asDe‘“*^  |+ 

83{-apApe‘“P*  -  a2A2e‘“2*  +  k^De*®*^  |  j 


(2.56) 

(2.57) 


(2.58) 


(2.59) 


(2.601 
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w=-i[a,{k,OpApe'“i'%k,G2A2e“‘J'+a,Ope‘“<"}+ 

a3{-apGpApe‘V_ajG2A2e“2'+k,GsDe‘“''}l 


Redefining  the  constants  such  that  the  coefficients  correspond  with  those  for  the 
displacement  field  in  the  upper  medium,  we  have 


-» Ap 


(2.62) 


-i  A2  =  *22 


2a 


(2.63) 


=  55  (2.64) 

where  <22,  and  %  are  the  displacement  potential  transmission  factors  for  the  Type  I, 
Type  n  and  shear  waves,  respectively.  The  displacements  in  the  porous  medium  can  then 
be  written 

+as'2'se‘“*^  }+ 

33  |-apTpe‘“p^  -  a  2T2e‘“2^ + k^T^e’®*^  |] 

^  J-l  (2.65) 

w(k„z)  =^[ai{k,TpGpe'“p%k^T2G2e‘“2^  }+ 

a3{-apTpGpe‘“p"  -  a2‘r2G2e‘“2"  +  k^T,Gsc‘“‘"  }] 


(2.66) 
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2.2.2  Wave  Velocities  and  Attenuations 

In  this  section  we  plot  Eq.  (2.46  and  2.52)  to  demonstrate  the  velocity  dependence 
of  the  Type  I,  Type  n,  and  shear  wave  on  frequency  and  sediment  type.  Three  sediment 
types  are  represented:  1)  an  unconsolidated  sand,  2)  an  unconsolidated  clay,  and  3)  a 
consolidated  sandstone.  The  material  properties  for  each  sediment,  required  by  Biot 
theory,  are  given  in  Table  2.1. 

Figures  2.3  and  2.4  show  the  frequency  dependence  of  the  phase  velocities  for  a 
marine  sand.  Material  properties  in  Table  2.1  and  were  derived  from  measurements  on  a 
shallow-water  carbonate  sand  [Badiey  et  al.  (1988)].  The  Type  I  wave  (the  classical  wave 
type  in  geophysics)  is  propagatory  at  all  frequencies  and  shows  a  phase  velocity  dispersion 
of  roughly  5%  over  this  frequency  range.  The  Type  n  wave  is  a  diffusive  wave  at  low 
frequencies  and  becomes  propagatory  at  nigh  frequencies  with  a  phase  velocity  about  one- 
tenth  of  that  for  the  Type  I  wave.  The  velocity  of  the  Type  n  wave  in  the  IdloHertz  region 
for  this  set  of  parameters  happens  to  be  close  to  that  measured  by  Paterson  (1956) 
(approximately  1 15  m/s)  on  a  fully  saturated  beach  sand  at  low  confining  pressures.  The 
shear  wave  velocity  predicted  by  Biot  theory  is  only  slightly  dispersive. 

Figures  2.5  -  2.7  show  the  attenuation  (scaled  by  frequency)  as  a  function  of 
frequency.  A  common  assumption  in  marine  sediment  acoustics  is  that  the  attenuation  is 
linearly  proportional  to  frequency  ( or  that  the  attenuation  scaled  by  frequency  is  constant). 
The  Biot  theory  predicts  a  non-linear  frequency  dependence  of  attenuation  due  to  the 
relative  motion  between  the  solid  matrix  and  the  interstitial  fluid.  At  low  frequency  the 
oscillating  fluid  obeys  Poiseuille  flow  where  the  velocity  distribution  across  a  (cylindrical) 
poTc.  is  parabolic,  with  zero  velocity  at  the  pore  wall  and  a  maximum  velocity  on  the  axis  of 
the  pore.  At  high  frequency,  the  fluid  obeys  Helmholtz  flow  where  the  fluid  inertia  causes 
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Table  2.1  Material  Properties  for  Three  Marine  Sediments 


Pronertv 

Units 

Sand 

Silty.-.Clay 

Sandstone 

Fluid  density,  (p^) 

g/cm^ 

1.024 

1.024 

1.024 

Ruid  bulk  modulus,  (Kj) 

dyne/cm^ 

2.38  IQlO 

2.38  lOlO 

2.38  lOlO 

Ruid  viscosity,  (il) 

poise 

.0101 

.0101 

.0101 

Grain  density,  (Pj) 

g/cm^ 

2.85 

2.68 

2.65 

Grain  bulk  modulus,  (K^) 

dyne/cm^ 

4.0  10»1 

3.5  10" 

3.7  10" 

Porosity,  (P) 

— 

.487 

.68 

.203 

Permeability,  (k) 

cm^ 

4.26  10-7 

5.2  10-10 

1.07  10-7 

Structiue  factor,  (c) 

— 

1.75 

3.0 

2.0 

Pore  size,  (a) 

cm 

4.18  10-3 

1.25  1(H 

2.5  10^ 

Frame  bulk  modulus 

(KJ 

dyne/cm^ 

3.33  108+ 

3.67  107+ 

/6.67  106 

8.32  1010+ 
jT.2  106 

Frame  shear  bulk  modulus 
(4) 

dyne/cm^ 

2.5  108+ 
/3.75  106 

7.86  106+ 
i2.5  105 

8.53  1010+ 
i4.25  108 
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the  velocity  distribution  to  be  nearly  flat  across  the  diameter  of  the  pore  with  an 
infinitesimal  boundary  layer  at  the  pore  wall.  The  frequency  separating  these  two  regimes 
is  termed  the  critical  frequency  and  is  defined  by 


r  _  _L 

57C  X  ■ 


(2.67) 


It  is  clear  from  Figures  2.5  -  2.7  that  only  far  below  or  above  the  critical  frequency 
(735  Hz)  is  the  attenuation  linearly  dependent  on  frequency  for  porous  media. 

Figures  2.8  through  2.14  show  the  frequency  dependence  of  velocities  and 
attenuations  for  a  marine  clay  and  sandstone.  The  influence  of  the  velocity  and  attenuation 
dispersion  shown  in  these  figures  on  seafloor  reflection  will  be  probed  in  detail  in  Chapter 
m.  Material  properties  for  these  sediments  are  given  in  Table  2.1;  the  silty  clay  properties 
were  taken  j5rom  Holland  and  Brunson  (1988)  and  the  properties  of  the  sandstone  were 
taken  from  Winkler  (1985). 


2.3  Application  of  the  Boundary  Conditions 

There  are  four  boundary  conditions  [Deresewiecz  and  Skalak  (1963)]  specified  at 
the  interface  between  a  fluid  and  a  poro-viscoelastic  solid: 

1 )  continuity  of  normal  fluid  displacement 

ftUo=Au-ft-w  ;z=0 

2)  continuity  of  normal  traction 

H(7u)-C(V.w)-2fA^  =  -p  ;  z  =  0  (2.68) 
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Figure  2.8  Type  I  wave  phase  velocity  for  silty  clay 


Figure  2.9  Type  11  and  shear  wave  phase  velocity  for  silty  clay 
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Figure  2.10  Type  I  wave  attenuation  for  silty  clay 
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Figure  2.14  Type  II  wave  attenuation  for  sandstone 
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equilibrium  of  tangential  traction 


;  z  =  0 


4)  equilibrium  of  fluid  pressure  (open  pores) 


M(V  •  w)  -  C(V  •  u  )=  p 


;  z  =  0 


and  the  Fourier  transforms  of  the  boundary  conditions  are: 


Uqz  ~ 


;  z  =  0 


3w, 


+2^ijk3iUx=-p 


(2.69) 


3ux  ..  . 

— ^-tkxUx,=0 

3)  ax  *  ^ 


Mr-ikxWx+^ 


]-c(^-ik,u,+^]=p 


Then,  substituting  the  equations  for  transformed  pressure  and  displacement  Eq  (2.17, 
2.22,  2.65,  2.66)  in  the  boundary  conditions,  the  expression  for  the  reflection  and 
transmission  factors  is  obtained  as: 
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a 

a„(l-G„) 

a^d  —  Go) 

Po 

CGp-H  2nk5 

CG2-H  .  2pk; 

0 

«P 

“2 

-Po 

C-MGp 

C-MGo 

w2 

v2 

1 

1 

O 

Vi 

'SR' 

a  ■ 

2k,a# 

o>^ 

-Po 

k 

•^2 

0 

0 

Po 

- 

(2.70) 


which  can  be  solved  (see  Appendix  B)  to  obtain  the  reflection  and  transmission  factors. 
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Chapter  3 

PLANE  WAVE  COEFFICIENTS  AND  SURFACE  WAVES 

3.1  Exact  Solution  for  Plane  Wave  Coefficients 

In  this  section,  the  plane  wave  reflection  coefficients  from  a  water- sediment 
halfspace  are  examined  where  the  sediment  is  modeled  as  a  poro-viscoelastic  solid,  as  a 
viscoelastic  solid,  and  as  a  fluid.  Both  consolidated  and  unconsolidated  sediments  are 
considered.  It  will  be  instructive  to  consider  the  reflectivity  from  simplest  model  of  the 
sediment  to  the  most  complex  (i.e.,  fluid,  solid,  then  poro-viscoelastic). 

3.1.1  Single  Phase  Media 

For  a  fluid-fluid  interface  the  reflectivity  is  governed  by 

-p  piVi  cos  9o  -  PoVqCOs  ei 
“  ”  plVi  cos  00  +  PoVo  cos  01  ^  ^ 


which  is  independent  of  frequency.  The  subscript  “o”  signifies  the  medium  where  the 
wave  originates  and  the  1  signifies  the  transmitting  medium.  For  a  non-lossy  material  the 
critical  angle  is  defined  as  the  angle  above  which  there  is  total  reflection  and  is  given  by 


;  vi>Vo 


(3.2) 


where  Vq  is  the  velocity  of  the  medium  where  the  wave  originates  and  Vi  is  the  velocity  of 
the  transmitting  medium.  For  any  real  material  with  frictional  and/or  viscous  losses,  there 
is  no  critical  angle,  since  there  is  no  total  reflection.  It  is,  however,  useful  to  refer  to  that 
angle  that  divides  angle  space  from  nearly  total  reflection  to  moderate  or  small  reflectivity. 
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This  angle  will  be  called  (for  lack  of  any  universally  applied  name)  the  pseudo-critical  (p-c) 
angle: 


.  .  ,  rRe{ki) 

0pc 


•]  ^ 


Re{ko}  >Re{ki} 


where  k  is  the  complex  wavenumber.  For  a  non-lossy  material  the  angle  of  intromission  is 
defined  as  that  angle  for  which  the  reflection  coefficient  is  zero  and  is  given  by: 


For  a  lossy  medium,  the  reflection  coefficient  is  non-zero  at  all  angles  and  the  pseudo¬ 
intromission  angle  is: 

1/2 

Re'lMZ  '.p'j  ]  :««(‘oi<R.(kl).PI>P„.  (3.5) 


Figure  3.1  shows  the  magnitude  of  the  plane  wave  reflection  coefficient  of  the 
amplitude  of  the  displacement  potential  for  a  sand,  silty  clay,  and  sandstone.  The  angle  of 
pseudo-intromission  exists  for  the  clay  and  pseudo-critical  angles  exist  for  the  sand  and 
sandstone.  The  geoacoustic  parameters  are  given  in  Table  3.1.  For  all  plots  in  this  chapter 
the  upper  halfspace  velocity  is  1524  m/s  and  the  density  is  1.024  g/cm^. 

For  a  fluid-viscoelastic  solid  interface  the  reflectivity  is  governed  by 


(3.6) 


piVi  cos  00 

cos2(2yi)  cos  Yi  +  ^sin2(2Yi)  cos  0i 

-  Po^o  COS01  cos  Yl 

piVi  cos  00 

cos2(2yi)  cos  Yl  +  ^sin2(2Yi)  cos  0i 

+  PoVo  COS01  cosYl 
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which  is  also  independent  of  frequency.  Vs  is  the  shear  velocity  and  Yl  is  the  angle  of  the 
refracted  shear  wave.  For  this  case  two  critical  angles  can  arise,  one  due  to  the 
compressional  wave  (as  in  the  fluid  case)  and  one  due  to  the  shear  wave; 


^cs 


=  sm 


-1 


Re  (ks)' 
Re  {ko)_ 


;Rc  {kol  >  Re  (ks). 


(3.7) 


Figure  3.2  shows  the  magnitude  of  plane  wave  reflection  coefficient  for  the  same 
marine  sand,  clay  ,  and  sandstone  case  as  in  Table  3.1  with  shear  wave  parameters  added 
(see  Table  3.2).  The  reflectivity  for  the  viscoelastic  sand  and  clay  sediments  is  essentially 
identical  with  that  for  die  fluid  sediments  (see  Figure  3.1)  indicating  that  the  conversion  of 
energy  to  shear  waves  is  negligible  for  these  material  parameters.  Loss  of  energy  to  shear 
waves  is  clearly  seen  in  the  sandstone  case  (compare  Figures  3.1  and  3.2)  and  the  shear 
pseudocritical  angle  is  apparent  at  5 1°.  The  geoacoustic  parameters  are  given  in  Table  3.2. 


Table  3.1  Geoacoustic  Parameters  for  the  Single-Phase  Fluid  Model 


Property 

Units 

Sand 

Silty-Clay 

Sandstone 

Compressional  velocity 

m/s 

1638.52 

1478.22 

3331.52 

Compressional  attenuation 

dB/trj/kHz 

.07 

.0018 

.039 

Density 

g/cm3 

1.961 

1.554 

2.320 

0  4 
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Figure  3.2  Plane  wave  reflection  coefficients  at  a  fluid-viscoelastic 
boundary 
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Table  3.2  Geoacoustic  Parameters  for  the  Single-Phase  Solid  Model 


Property 

Comprcssional  velocity 
Compressional  attenuation 
Shear  velocity 
Shear  attenuation 
Density 


Units 

Sand 

m/s 

1638.52 

dB/m/kHz 

.070 

m/s 

121.83 

dB/m/kHz 

4.474 

g/cm3 

1.961 

Silty-Clay 

Sandstone 

1478.22 

3331.52 

.0018 

.039 

22.5 

1949.65 

38.83 

.239 

1.554 

2.320 

3.1.2  Two-Phase  Porouc  Media 

For  a  fluid-poro-viscoelastic  interface,  the  plane  wave  reflection  and  transmission 
coefficients  are  related  to  the  reflection  and  transmission  factors  in  Eq.  (2.70)  by 

(3-8) 

Tp  =  ‘2pe»«h 

=  *22  e‘ 

Ts  =  'rse»ah 


Substituting  these  definitions  in  Eq.  (2.70),  the  factor  drops  out  and  the  plane  wave 
coefficients  depend  solely  on  material  properties  and  not  on  source  height 
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3.1.2.1  Sand:  Vp  >  Vo,  Vj  <  Vo 

Figure  3.3  shows  the  plane  wave  reflection  coefficient  of  the  amplitude  of  the 
displacement  potential  for  the  marine  sand  given  in  Table  2.1.  There  are  three  important 
differences  between  the  reflectivity  from  porous  material  and  reflectivity  from  a  single 
phase  fluid  or  solid  material.  First,  note  that  reflectivity  from  a  poro-viscoelastic  medium 
is  frequency  dependent  In  particular  the  p-c  angle  is  a  function  of  frequency  and  decreases 
monotonically  with  increasing  frequency  (from  roughly  82°  at  10  Hz  to  70°  at  10  kHz).  It 
is  clear  that  the  decrease  in  the  p-c  angle  arises  from  the  increase  in  the  phase  velocity  of  the 
Type  I  wave  (which  is  the  only  wave  type  that  meets  the  criterion  Re{ko}  >  Re  (ki)). 
Figure  3.4  shows  the  frequency  dependence  of  the  p-c  angle  based  on  Eq.  (3.2)  for  the 
marine  sand  parameters  in  Table  2.1. 

The  second  feature  of  importance  in  Figure  3.3  is  the  dependence  of  the  reflection 
coefficient  on  frequency  above  the  critical  angle.  In  this  angular  regime,  the  reflection 
coefficient  is  inversely  proportional  to  the  attenuation  per  cycle.  The  minimum  of  the 
reflection  coefficient  occurs  at  the  peak  of  the  attenuation  (see  Figure  2.5)  between  100  Hz 
and  1000  Hz. 

The  third  feature  of  interest  is  the  reflectivity  at  normal  incidence.  For  a  single 
phase  fluid,  the  reflection  coefficient  at  ntnmal  incidence  is 


Rff(0o=O)  = 


PlVl  -  PqVq 
piVi  +poVo 


(3.9) 
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Figure  3.3  Plane  wave  reflection  cwfficients  at  a  fluid-poro-viscoelastic 
boundary.  The  sediment  is  sand. 


Pa«udo-CriUc«l  AngU  (dag) 
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which  increases  as  Vi  increases.  Since  we  have  seen  (e.g.  Figure  2.3)  that  for  porous 
media  Vi  increases  as  a  function  of  frequency  we  might  guess  that  normal  incidence 
reflectivity  increases  as  a  function  of  frequency.  From  Figure  3.3,  however,  reflectivity 
decreases  as  a  function  of  ticquency.  This  is  due  to  the  excitation  of  Type  n  waves  in  the 
porous  medium  as  frequency  increases  as  shown  in  Figure  3.5. 

3.1.2.2  Silty  Clay:  Vp  <  Vo,  Vg  <  Vo  ' 

Figure  3.6  shows  the  reflectivity  from  a  fluid-silty  clay  interface  where  the 
parameters  for  the  silty  clay  are  given  in  Table  2.1.  Over  this  frequency  band  the  sediment 
has  a  negligible  velocity  dispersion  (see  Figure  2.8)  and  a  nearly  linear  dependence  of 
attenuation  (see  Figure  2.10)  on  frequency,  thus  the  reflectivity  is  also  nearly  independent 
of  frequency.  Also,  the  shear  wave  and  Type  n  wave  in  the  clay  have  a  low  velocity  with 
respect  to  that  of  the  incoming  wave  so  that  the  reflectivity  is  negligibly  influenced  by  these 
wave  types. 

3. 1.2.3  Sandstone:  Vp  >  Vo,  Vj  >  Vo 

Figure  3.7  shows  the  reflection  coefficient  for  the  sandstone  parameters  given  in 
Table  2.1.  The  critical  angle  for  the  compressional  wave  (26°)  and  the  shear  wave  (50°) 
are  both  apparent.  Reflectivity  decreases  with  frequency  monotonically  at  all  angles  of 
incidence  due  to  the  increasing  transmissivity  of  the  Type  n  wave  (see  Figure  3.8).  The  p- 
c  angle  for  the  Type  I  wave  is  nearly  independent  of  frequency  since  the  ratio  Vp/Vo 
changes  less  than  .2%  over  this  frequency  band.  The  p-c  angle  for  the  shear  wave, 
however,  appears  to  decrease  with  increasing  frequency.  The  shear  velocity  is  indeed 
dispersive,  but  increases  with  frequency  so  that  from  Eq.  (3.7)  we  would  expect  the  p-c 
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Plane  wave  reflection  coefflcients  at  a  fluid-poro-viscoelastic 
boundary.  Tbe  sediment  is  silty-clay. 
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Figure  3.8  Plane  wave  transmission  coefficients  for  the  Type  n  wave  at  a 
fluid-poro-viscoelastic  boundary.  The  sediment  is  sandstone. 
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angle  to  increase  slightly  with  angle.  The  decrease  in  reflectivity  around  the  shear  p-c  angle 
(which  manifests  itself  as  an  apparent  shift  in  p-c  angle)  is  due  instead  to  the  increase  of 
transmissivity  of  the  Type  n  wave  near  this  angle  (see  Hgure  3.8). 

3.2  Implications  for  Acoustic  Propagation  Predictions 

In  support  of  acoustic  propagation  experiments,  sediment  cores  are  commonly 
taken  to  obtain  measurements  of  compressional  velocity  (occasionally  also  shear  velocity) 
and  attenuation.  Recently,  techniques  have  been  developed  to  measure  these  quantities  in 
situ.  Typical  frequencies  for  either  lab  or  in  situ  measurements  are  in  the  hundreds  of 
kilohertz.  In  developing  a  low  frequency  geoacoustic  model,  the  velocity  is  generally 
assumed  independent  of  frequency  and  the  attenuation  is  assumed  linearly  dependent  on 
frequency  and  extrapolated  down  to  the  frequency  of  interest  In  this  section  the  effects  of 
this  oversimplified  extrapolation  arc  demonstrated,  by  comparing  reflectivity  predicted  by 
current  methodology,  as  outlined  above,  and  Biot  theory. 

3.2.1  Sand;  Vp  >  Vo,  Vj  <  Vq 

Figure  3.9  shows  the  reflectivity  from  a  poro-viscoelastic  solid  as  a  function  of 
frequency  along  with  the  prediction  of  a  fluid  geoacoustic  model  that  has  been  extrapolated. 
(The  velocity  and  attenuation  of  the  fluid  approximation  to  the  sediment  were  calculated 
from  Biot  theory  at  400  kHz  where  measurements  of  these  properties  are  most  commonly 
taken).  The  differences  between  the  fluid  sediment  model  and  porous  sediment  model  are 
striking.  For  long-range  propagation  in  shallow  water  the  reflectivity  at  high  angles  tends 
to  determine  the  characteristics  of  the  transmission  loss:  the  location  of  the  p-c  angle 
governs  the  number  of  modes  that  will  propagate  and  the  reflectivity  above  the  p-c  angle 
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governs  the  modal  attenuation  coefficient  The  extrapolated  fluid  model  thus  overpredicts 
the  number  of  propagating  modes  and  underpredicts  the  modal  attenuation  coefiicient  for 
propagation  between  100  Hz  and  several  kHz.  The  results  for  extrapolating  a  solid  model 
from  400  kHz  would  be  essentially  identical  to  that  of  the  fluid  shown  in  Figure  3.9. 

3.2.2  Silty  Clay:  Vp  <  Vo,  Vs  <  Vo 

Figure  3.10  shows  the  reflectivity  fiom  a  poro- viscoelastic  silty  clay  as  a  function 
of  frequency  along  with  the  prediction  of  the  extrapolated  fluid  model  (predictions  of  the 
extrapolated  fluid  and  exO'apolated  solid  model  were  identical  to  within  the  width  of  a 
plotted  line).  The  complex  velocity  for  the  extrapolated  model  was  obtained  from  Biot 
theory  at  4(X)  kHz  as  in  the  sand  case.  The  predictions  of  reflectivity  using  the  single-phase 
fluid  approximation  are  quite  comparable  to  the  full  Biot  predictions  except  near  the  angle 
of  intromission,  where  the  full  theory  predicts  a  lower  reflectivity.  At  the  angle  of 
intromission,  the  reflection  coefGcient  is  proportional  to  the  compressional  attenuation  per 
wavelength;  Biot  theory  predicts  a  non-linear  fiequency  dependence  on  attenuation  between 
400  kHz  and  the  kilohertz  regime  (see  Figure  2.10).  Thus,  the  attenuation  for  the 
extrapolated  model  is  too  high. 

Seafloor  reflection  experiments  confirm  that  measurements  of  compressional  wave 
attenuation  at  400  kHz  caiuiot  be  linearly  extrapolated  down  to  the  50  Hz  -  1500  Hz  band. 
In  this  band  the  strongest  arrival  from  an  explosive  source  impinging  on  the  seafloor  at  low 
grazing  angles  comes  from  paths  that  refract  through  the  sediment  accumulating  a  loss 
proportional  to  sediment  attenuation.  Estimates  of  effective  attenuation  from  these  paths 
range  fi'om  .001  dB/m/kHz  -  .03  dB/m/kHz  (including  measurements  at  roughly  200  deep- 
ocean  sites  in  the  North  Atlantic,  Pacific,  Mediterranean  and  Indian  Ocean  (C.W.  Spofford, 
personal  communication)).  Measured  intrinsic  attenuation  in  the  hundreds  of  kHz  across  a 
core  barrel  for  deep  sea  sediment  types,  however,  has  been  in  the 
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Figure  3.9  Comparison  of  reflectivity  from  a  commonly  used  geoacoustic 
extrapolation  technique  with  Biot  theory.  The  sediment  is 
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range  .03  dB/m/kHz  -  .10  dii/m/kHz.  (see  Hamilton  (1980)).  Since  effective  attenuation 
must  be  greater  than  or  equal  to  intrinsic  attenuation  it  is  clear  that  the  experimental  evidence 
points  to  a  non-linear  frequency  dependence  of  attenuation  over  the  band  50  Hz  -  400  kHz. 
3.2.3  Sandstone:  Vp  >  Vq,  Vs  >  Vq 

Figure  3.1 1  shows  the  reflectivity  from  a  poro- viscoelastic  sandstone  as  a  function 
of  frequency  along  with  the  prediction  of  the  extrapolated  solid  and  extrapolated  fluid 
model.  Again,  the  complex  velocities  for  the  extrapolated  models  were  obtained  at  400 
kHz.  As  expected,  the  fluid  model  gives  a  very  poor  approximation  to  the  poro- 
viscoelastic  results  since  no  shear  wave  p-c  angle  is  predicted.  The  extrapolated  solid 
model  predistions  are  tolerable  at  the  low-frequencies  (less  than  10  kHz)  but  at  high 
frequencies  the  transmissivity  of  the  Type  n  wave  can  reduce  reflectivity  by  as  much  as 
50%. 

3.3  Approximate  Models 

Several  attempts  have  been  made  to  use  Biot  theory  in  a  simplified  “Biot  fluid” 
formulation  for  unconsolidated  sediments.  This  should  be  a  reasonable  approximation  at 
low  frequencies  since  tire  Type  n  wave  has  a  negligible  influence  (at  least  at  a  single 
boundary)  and  the  role  of  the  sediment  rigidity  also  is  negligible.  One  such  reflection 
model  might  be  Eq.  (3.1)  with  a  frequency-dependent,  complex  velocity  obtained  from 
Biot  predictions  for  the  Type  I  wave.  Likewise  a  "Biot  solid"  approximation  model  might 
be  Eq  (3.6)  with  the  velocities  obtained  from  Biot  predictions  for  the  Type  I  wave  and  the 
shear  wave.  Figure  3.12  contains  such  predictions  compared  with  the  full  theory.  In 
Figtire  3.12,  the  "Bin?  fluid"  approximation  to  the  carbonate  sand  is  quite  good  at  low 
frequencies  for  all  angles  of  incidence.  The  reflectivity  under  the  approximate  model, 
however,  increases  with  increasing  frequency  near  and  at  normal  incidence;  the  full  theory 
predicts  a  decreasing  reflectivity  with  increasing  frequency  due  to  increased  coupling  to  the 
Type  n  wave.  At  low  frequencies  (below  100  Hz)  the  results  of  the  full  and  approximate 
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3.11  Comparison  of  reflectivity  from  a  commonly  used  geoacoustic 
extrapolation  technique  with  Biot  theory.  The  sediment  is 
sandstone. 
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Figure  3.12a,b 


Comparison  of  Biot  theory  with  an  approximate  model  (#). 
The  sediment  is  sand,  the  frequency  is:  a)  10  Hz;  b)  100 


Figure  3.12c, d  Comparison  of  Biot  theory  with  an  approximate  model  (#). 

The  sediment  is  sand,  the  frequency  is:  c)  1  kHz;  d)  10 
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model  are  essentially  identical  near  normal  incidence.  At  10  kHz,  the  approximate  model 
predicts  reflectivity  about  15%  Iiigher  than  dic  full  theory  near  normal  incidence.  In  Figure 
3.13,  the  "Biot  solid"  approximation  to  the  sandstone  is  quite  good  at  low  frequencies 
(below  1  kHz).  At  higher  frequencies,  the  Type  II  wave  goes  from  diffusive  to 
propagatory  and  more  energy  is  transmitted  to  this  wave  type.  Therefore,  the  approximate 
model  that  does  not  include  this  wave  type  is  no  longer  applicable,  and  the  full  theory  is 
required. 

Beebe  (1981)  proposed  an  even  simpler  model,  where  the  “Biot  fluid”  is  modeled 
as  having  a  constant  phase  velocity  and  an  attenuation  where  the  frequency  exponent  is 
calculated  by  Biot  theory  over  the  band  of  interest  (in  this  case  50  Hz  -  1000  Hz).  The 
U.S.  Navy  standard  low  frequency  bottom  loss  data  base  later  employed  a  similar  sediment 
model,  where  the  frequency  exponent  is  obtained  by  inversion  of  TL  data  rather  than  Biot 
theory.  It  is  crucial,  in  this  approximate  model,  that  the  constant  phase  velocity  is  obtained 
in  the  frequency  band  of  interest.  For  example,  in  Beebe's  study  the  velocities  were 
obtained  by  dispersion  analysis  in  the  same  frequency  band  in  which  the  mode  attenuation 
coefficients  were  measured  resulting  in  ideal  conditions  for  the  approximate  model  to 
produce  good  results.  A  similar  condition  is  obtained  synthetically  in  Figure  3.14  by  using 
a  constant  velocity  (as  predicted  by  Biot  theory  in  the  middle  of  the  band  of  interest  (50  Hz 
-  1000  Hz).  Figure  3.14  compares  the  approximate  model  with  the  full  theory  (for  the 
sand)  where,  in  the  approximate  model,  the  imaginary  part  of  the  phase  velocity  comes 
from  Biot  theory  and  the  real  part  is  a  constant  value  commensurate  with  the  center  of  the 
band.  The  results  are  quite  reasonable  although  the  approximate  model  breaks  down  with 
increasing  frequency  as  the  assumed  velocity  becomes  a  poorer  approximation  of  the  true 
velocity. 

Figure  3.15  shows  the  same  comparison  if  the  approximate  model  were  developed 
using  velocity  data  from  core  analysis.  That  is,  the  velocity  used  in  the  approximate 
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Figure  3.13  Comparison  of  Biot  theory  with  an  approximate  model  (#) 
The  sediment  is  sandstone. 


Angle  of  Incidence  (deg) 


Figure  3.14  Comparison  of  Biot  theory  with  the  constant  phase  velocity 
approximate  model  (*).  The  constant  phase  velocity  was 
derived  within  the  50  Hz  •  1000  Hz  band.  Frequency  is:  a)  SO 
Hz;  b)  1000  Hz. 
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Figure  3.15  Comparison  of  Biot  theory  with  the  constant  phase  velocity 
approximate  model  (*).  The  constant  phase  velocity  is  derived 
at  400  kHz  commensurate  with  velocity  measurements  on 
sediment  cores.  Frequency  is:  a)  50  Hz;  b)  1000  Hz.. 
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model  here  was  obtained  at  400  kHz  from  Biot  theory.  It  is  clear  that  using  core  data  for 
input  to  the  approximate  model  yields  very  poor  results. 

3 . 4  Interface  Waves 


Before  proceeding  to  the  inverse  transform  we  consider  the  interface  waves  at  a 
fluid-poro-viscoelastic  boundary,  the  properties  of  which  may  be  ascertained  from  the  same 
set  of  simultaneous  equations  that  govern  the  plane  wave  reflection  and  transmission 
coefficients.  Various  researchers  have  investigated  related  surface  waves;  however,  it 
appears  that  the  general  case  of  surface  waves  at  a  fluid-poro-viscoelastic  boundary  (with 
complex  wave  speeds)  has  not  been  addressed.  Deresewicz  (1964)  derived  the  secular 
equation  without  accounting  for  losses  due  to  friction  at  the  grain-to-grain  contacts  and 
Feng  and  Johnson  (1983a  and  1983b)  considered  the  high  frequency  limit  to  the  Biot 
theory  where  the  viscous  attenuation  is  zero. 

The  dispersion  relation  for  these  waves  can  be  obtained  by  setting  the  denominator 
of  the  plane  wave  reflection  and  transmission  coefficients  to  zero.  From  Cramers  Rule, 
and  Eq.  (2.70)  we  have  that 
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where  the  surface  wave  velocity,  Vgurf  =  Three  kinds  of  surface  waves  can  occur  in 
the  solution;  a  true  surface  wave  that  travels  slower  than  all  the  wave  speeds,  a  pseudo- 
Scholte  wave  that  travels  between  the  shear  wave  velocity  and  the  Type  11  wave  velocity 
(that  leaks  energy  to  the  slower  Type  II  wave),  and  a  pseudo-Rayleigh  wave  that  occurs 
when  the  surface  wave  velocity  exceeds  the  velocity  of  the  overlying  water  column. 

Mathematically,  the  pseudo-Rayleigh  wave  arises  from  a  pole  on  the  lower 
Reimann  sheet  where  Im{o}  >  0  and  Iin{a2}  >  0 .  The  tme  surface  wave  is  located  on  the 
top  Reimann  sheet  where  Im{a  J  <  0,  and  the  amplitude  decays  exponentially  away  from 

the  boundary.  The  pseudo-Rayleigh  wave  amplitude  actually  increases  away  from  the 
boundary  (in  the  fluid)  and  this  occurs  because  the  wave  radiates  into  the  fluid  at  its  trace 
velocity  so  that  in  a  given  plane  perpendicular  to  the  boundary,  the  points  closer  to  the 
boundary  in  the  fluid  receive  contributions  from  the  wave  further  along  the  boundary 
(which  are  less  energetic).  TTiis  behavior  of  the  received  amplitude  of  the  pseudo-Rayleigh 
wave  is  entirely  commensurate  with  the  behavior  of  a  lateral  wave. 

The  complex  surface  wave  velocity  from  a  single-phase  solid  will  be  calculated  in 
order  to  compare  with  the  results  from  poro-viscoelastic  theory.  The  secular  equation  is 
given  by  Brekovskikh  (1980)  as: 
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3.4.1  Sand:  Vp  >  V©,  Vs  <  V® 

For  the  same  sand  case  considered  previously,  the  surface  wave  velocity  as  a 
function  of  frequency  is  plotted  in  Figure  3.16  along  with  the  phase  velocity  of  the  Type  II 
and  shear  wave.  Note  that  the  surface  wave  changes  continuously  from  a  pseudo-Scholte 
wave  to  a  true  surface  wave  at  about  700  Hz^.  Also  plotted  in  Figure  3.16  is  the  surface 
wave  velocity  in  the  single  phase  solid  approximation,  where  the  complex  phase  velocities 
used  to  calculate  the  surface  wave  velocity  were  taken  from  Biot  theory  at  400  kHz.  The 
single  phase  approximation  to  the  surface  wave  velocity  is  roughly  10%  higher  than  the  full 
poro-viscoelastic  predictions  over  this  frequency  band.  The  attenuation  of  the  poro- 
viscoelastic  surface  wave  shows  a  similar  frequency  dependence  to  the  Type  I  and  shear 
wave  (Figure  3. 17).  The  single-phase  solid  approximation  to  the  surface  wave  attenuation 
is  in  general  too  conservative. 

3.4.2  Silty  Clay:  Vp  <  V®,  Vj  <  Vp 

The  surface  wave  velocity  for  the  silty  clay  changes  continuously  from  a  pseudo- 
Scholte  wave  to  a  true  surface  wave  at  167  kHz  (see  Figure  3.18).  The  velocity  dispersion 
of  the  surface  wave  is  small  so  that  the  surface  wave  velocity  under  ihc  single-phase  solid 
approximation  is  quite  good  for  this  case.  The  single-phase  approximation  to  the 
attenuation  (Figure  3.19),  however,  is  very  poor. 


1.  The  surface  wave  velocity  and  attenuation  at  4.5  Hz  (101  m/s  and  25  dB/km)  is 
coincidentally  comparable  to  that  measured  by  Shirmer  (1980)  (111  rn/s  and  7  dB/km)  on  a 
sandy  bottom  (of  unknown  material  parameters)  in  the  North  Sea. 
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Fieure  3.16  Phase  velocity  of  various  wave  types  in  carbonate  sand, 
pseado-Scholte  wave  at  a  fluid-poro-viscoelastic  boundary 
(porous  SW);  pseudo-SchoIte  wave  at  a  fluid-viscoelastic 
boundary  (solid  SW);  shear  wave  in  a  poro-viscoelastic 
sediment  (shear);  and  Type  II  wave  in  a  poro-viscoelastic 
sediment  (Type  II). 
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Figure  3.17  Attenuation  of  various  wave  types  in  carbonate  sand;  pseudo* 
Scholte  wave  at  a  fluid<poro-viscoeIastic  boundary  (porous 
SW);  pseudo-Scbolte  wave  at  a  fluid*viscoelastic  boundary 
(solid  SW);  and  shear  wave  in  a  poro* viscoelastic  sediment 
(shear). 
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Figure  3.18  Phase  velocity  of  various  wave  types  in  silty  clay;  pseudo* 
Scbolte  wave  at  a  fluid-poro-viscoelastic  boundary  (porous 
SW);  pseudO'Scholte  wave  at  a  fluid-viscoelastic  boundary 
(solid  SW);  shear  wave  in  a  poro-viscoelastic  sediment 
(shear);  and  Type  II  wave  in  a  poro-viscoelastic  sediment 
(Type  II). 
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Figure  3.19  Attenuation  of  various  wave  types  in  silty  clay;  pseudo-Scholte 
wave  at  a  fluid-poro- viscoelastic  boundary  (porous  SW); 
pseudo-Scbolte  wave  at  a  fluid-viscoelastic  boundary  (solid 
SW);  and  shear  ^ave  in  a  poro-viscoelastic  sediment  (shear). 
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3.4.3  Sandstone:  Vp>Vo,  Vs>Vo 

For  this  sediment  (Vs>Vo),  both  the  pseudo-SchoIte  wave  and  the  pseudo- 
Rayleigh  wave  exist.  For  both  the  pseudo-Scholte  (Figures  3.20  and  3.21)  and  the 
pseudo-Rayleigh  (Figures  3.22  and  3.23)  surface  wave  types,  the  single-phase 
approximation  is  reasonable  for  the  velocity  but  quite  poor  for  the  attenuation. 

There  is  some  difficulty  associated  with  assigning  a  wave  speed  and  attenuation  to 
the  pseudo-Rayleigh  wave.  Consider  two  studies  of  a  fluid-solid  boundary  with  surface 
waves  excited  from  an  impulsive  source.  Van  der  Hijden  (1984)  states  that  “from  a 
mathematical  point  of  view  it  is  not  appropriate  to  associate  with  [Vsurf]  a  pseudo-Raylei£.i 
wave  speed  because  there  is  no  wave  front  associated  with  this  phenomenon.”  He 
considers  the  possibility'  of  using  the  true  Rayleigh  wave  speed  (calculated  as  if  a  vacuum 
replaced  the  fluid)  to  describe  the  velocity  of  a  pulse  and  attributing  the  change  in  shape  of 
the  pulse  to  some  dispersive  property.  He  opts  instead  to  use  Vsurf  as  the  “pseudo- 
Rayleigh  travel  speed”  which  describes  the  speed  of  the  pulse  maximum.  Roever  et  al. 
(1959)  on  the  other  hand  demonstrate  experimentally  and  theoretically  that  the  zero¬ 
crossing  of  a  pseudo-Rayleigh  wave  does  travel  along  the  boundary  at  precisely  the 
Rayleigh  wave  velocity;  that  is  independent  of  the  fluid  properties!  Roever  et  al.  (1959) 
indicate,  however,  that  the  zero  crossing  of  a  pseudo-Rayleigh  wave  excited  by  a  time 
harmortic  source  does  travel  at  Vsurf  (i.e.,  is  influenced  by  the  properties  of  the  upper 
halfspace).  Since  this  study  considers  the  response  of  the  boundary  of  a  time  harmonic 
source.  Figures  3.22  and  3.23  represent  solutions  to  the  secular  equation  accounting  for  the 
properties  of  the  upper  halfspace. 
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Figure  3.20  Phase  velocity  of  various  wave  types  in  sandstone.  Pseudo- 
Scholte  wave  at  a  fluid-poro-viscoelastic  boundary  (porous 
SW);  pseudo-SchoIte  wave  at  a  fluid-viscoelastic  boundary 
(solid  SW);  and  shear  wave  in  a  poro-viscoelastic  sediment 
(shear). 
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Figure  3.21  Attenuation  of  various  wave  types  in  sandstone;  pseudo- 
Scholte  wave  at  a  fluid-poro<viscoelastic  boundary  (porous 
SW);  pseudo-Scholte  wave  at  a  fIuid>viscoelastic  boundary 
(solid  SW);  and  shear  wave  in  a  poro-viscoelastic  sediment 
(shear). 
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Figure  3.22  Phase  velocity  of  various  wave  types  in  sandstone.  Pseudo* 
Rayleigh  wave  at  a  fluid-poro-viscoelastic  boundary  (porous 
RW);  pseudo-Rayleigh  wave  at  a  fluid* viscoelastic  boundary 
(solid  RW);  and  shear  wave  in  a  poro*viscoelastic  sediment 
(shear). 
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Figure  3.23  Attenuation  of  various  wave  types  in  sandstone;  pseudo* 
Rayleigh  wave  at  a  fluid-poro*viscoelastic  boundary  (porous 
RW);  pseudO’Rayleigh  wave  at  a  fluid-viscoelastic  boundary 
(solid  RW);  and  shear  wave  in  a  poro-viscoelastic  sediment 
(shear). 
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3.5  Summary  and  Conclusions 

In  this  chapter,  plane  wave  coefficients  for  fluid-fluid,  fluid-viscoelastic,  and  fluid- 
poro-viscoelastic  boundaries  were  compared.  Biot  theory  predicts  frequency  dependent 
plane  wave  reflection  (and  transmission)  coefficients  as  compared  with  the  frequency 
independent  coefficients  predicted  using  single-phase  fluid  or  viscoelastic  theories. 

For  the  sand  case  the  critical  frequency  (Eq.  2.67)  was  low  enough  to  be  in  the 
ftcquency  range  of  interest  Therefore,  the  velocity  and  attenuation  of  each  wave  type 
(especially  the  Type  I  wave  which  is  the  most  important  wave  influencing  reflectivity  for 
this  sediment  type)  varied  non-linearly  over  the  band  of  interest  leading  to  the  frequency 
dependent  plane  wave  coefficients.  For  the  silty  clay  case  the  critical  frequency  was  high 
enough  that  the  velocity  of  the  Type  I  wave  (the  only  wave  participating  strongly  in  the 
reflection  and  transmission  process)  was  nearly  constant  with  fiiequency  and  the  attenuation 
was  slightly  non-linear  with  frequency.  The  reflection  coefficient  was  independent  of 
frequency  over  the  10  Hz  -  10  kHz  band.  For  the  sandstone  case,  the  velocities  of  the 
Type  I  and  shear  waves  were  only  weakly  dispersive.  From  previous  results  this  suggests 
that  the  plane  wave  coefficients  should  exhibit  only  a  weak  frequency  dependence.  In  this 
case,  however,  the  velocity  dispersion  of  the  Type  n  significantly  influenced  the  plane 
wave  coefficients  at  high  fiequency. 

In  attempting  then  to  specify  criteria  for  when  the  single-phase  models  may  be 
substituted  for  the  full  theory,  two  factors  must  be  taken  into  account:  1)  the  frequency 
dependence  of  the  plane  wave  coefficients  (arising  from  the  frequency  dependence  of  the 
complex  wave  velocities)  and  2)  the  role  of  the  Type  II  waves. 

First  consider  the  frequency  dependence.  It  is  the  critical  frequency  that  determines 
whether  or  not  the  full  Biot  theory  is  required  for  predicting  energy  balance  at  the  seafloor. 
If  the  critical  frequency  for  a  particular  sediment  type  is  well  outside  the  frequency  band  of 
interest,  then  a  single  phase  model  will  give  a  good  approximation  to  the  reflection 
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coefficient.  If  the  critical  frequency  (Eq.  (2.67))  is  near  the  frequency  band  of  interest, 
however,  the  complex  wave  velocities  will  vary  substantially  across  that  band  and  the  full 
Biot  theory  will  be  required.  This  last  conclusion  presupposes  that  the  reflection  equation 
will  be  sensitive  to  the  frequency  dependence  of  at  least  one  of  the  three  wave  speeds 
predicted  by  Biot  theory.  While  this  may  not  be  true  for  an  arbitrary  fluid-  porous  material 
interface,  it  is  expected  to  be  true  for  the  ranges  of  interface  types  found  at  the  seafloor.  A 
word  is  in  order  here  regarding  the  phrase:  “frequency  band  of  interest”.  Certainly  it  refers 
to  the  frequencies  over  which  the  acoustic  measurements  were  taken.  However,  if 
geoacoustic  properties  were  measured  on  core  samples  at  high  frequencies,  the  “frequency 
band  of  interest”  referred  to  in  the  above  criterion  must  be  expanded  to  include  these 
frequencies.  It  was  shown  that  the  traditional  extrapolation  of  geoacoustic  properties  from 
core  data  can  lead  to  erroneous  results  (at  least  with  respect  to  Biot  theory)  when  the  critical 
frequency  is  between  the  measurement  bands  of  the  acoustic  data  and  the  geoacoustic  data. 

The  second  factor  the  role  of  the  Type  11  wave  is  somewhat  simpler.  The  Type  II 
waves  are  only  important  to  reflection/transmission  at  high  frequencies  where  they  arc 
propagatory:  “high  frequencies”  roughly  corresponds  to  frequencies  above  the  critical 
frequency. 

The  surface  waves  arising  at  a  fluid-poro-viscoelastic  boundary  were  also  examined 
and  compared  with  the  classical  single-phase  solid  approximation.  The  single  phase  solid 
approximation  provided  a  reasonable  approximation  to  the  surface  wave  phase  velocity  (as 
predicted  by  Biot  theory)  for  all  three  sediment  types,  however,  in  all  cases  neither  the 
magnitude  nor  the  frequency  dependence  of  the  attenuation  was  adequately  predicted. 
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Chapter  4 

THEORY  n  -  INVERSE  TRANSFORM 


In  this  chapter  the  integral  expressions  for  the  field  quantities  in  the  upper  and  lower 
halfspace  are  evaluated  by  a  “mixed”  approach,  a  numerical  integration  and  by  the  Fast 
Fourier  Transform  (FFT).  The  mixed  approach  consists  of  breaking  the  field  into 
components  representing,  for  example,  contributions  from  the  direct  wave,  the 
reflected/transmitted  geometrical  optics  wave,  and  the  various  head  and  surface  waves. 
Results  from  the  mixed  approach  and  the  FFT  are  compared  to  examine  the  bounds  of 
applicability  for  each. 

4.1  Mixed  Approach 

4.1.1.  Field  in  Upper  Halfspace 

The  expression  for  the  pressure  field  in  the  upper  halfspace  is  obtained  by  Fourier 
Transforming  Eq.  (2.17): 


p(x.z2  0)=-i^ 
47C 


7  {e“‘“«^-‘’l  +  9le-**"}a-*e"‘^^dk^ 


(4.1) 


where  the  first  term  corresponds  to  the  direct  field  and  the  second  term  to  the  reflected  field. 
The  reflected  field  can  be  converuently  decomposed  into  (1)  a  contribution  from  plane 
waves  reflected  at  the  interface  near  angles  predicted  by  Snell's  law  (the  geometrical  ray 
term);  (2)  a  contribution  from  head  waves  (arising  from  branch  points  in  the  reflection 
factor  91);  and  (3)  a  contribution  from  surface  waves  (arising  from  poles  in  9i) .  Each  of 
these  contributions  will  be  considered  separately  in  the  follow;  ng  sections. 
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4. 1.1.1  The  Direct  Field 

The  integral  expression  for  the  direct  field  can  be  solved  exactly: 


P(l(x,2  >  0)  =  +  (z  -  j 


which  represents  an  outgoing  cylindrical  wave  radiated  fiom  the  line  source. 


4. 1.1.2  Reflected  Field  (Geometrical  Ray  Path) 
The  integral  expression  for  the  reflected  field  is 


where  the  plane  wave  reflection  coefficient  has  been  substituted  for  the  reflection  factor 
using  Eq.  (3.8).  Making  the  transformation,  kx  =  ko  sin  0,  and  employing  the  coordinate 
system  (see  Figure  4.1):  x  =  s  sin  Gq  and  z+h  =  s  cos  0o,  Eq.  (4.3)  can  be  written: 


Pr(x,z>0)  = 


2  2 


The  method  of  steepest  descent  is  employed  to  evaluate  this  integral.  The  method 


considers  integrals  of  the  form 


I=LF(0)eP^(®^d0 


Figure  4.1  Coordinate  system  for  the  reflected  field, ’h*  is  the  source 
height. 
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in  the  complex  plane,  where  F(0)  must  be  a  slowly  varying  function  of  0.  The  technique 
consists  of  finding  the  saddle  point  of  the  function  f(0)  and  the  path  along  which  the 
function  decreases  most  rapidly.  The  original  path  of  integration  is  then  deformed  along 
that  path  and  the  integrand  is  expanded  in  a  Taylor  series  around  the  saddle  point.  Since 
most  of  the  contribution  to  the  integral  comes  from  near  the  saddle  point,  the  function  can 
be  approximated  by  a  few  terms  and  the  integral  is  evaluated  term  by  term. 

Letting  F(0)  =^,  p  =  Iq,  s,  and  f(0)  =  -t  cos  (0-0o),  the  saddle  point,  0sp,  is  defined  at 
f(0)  =  0  which  gives 


0sp  — 


(4.6) 


Thus,  the  saddle  point  is  simply  located  at  the  angle  of  incidence  corresponding  to  the 
geometrical  ray  path  and  is  real  for  real  velocity  Vq,  varying  from  0  to  ?r/?  radians.  T^e 
path  of  steepest  descent  is  defined  where  the  imaginary  pan  of  f(0)  is  constant. 


Im{f(9)}=Im{f(9^)} 


(4.7) 


or 


Re{cos0cos0o  +sin0sin0o}=l 


Letting  0=0^  +  i0l  the  path  of  steepest  descent  is  defined  by. 


(4.8) 


cosh0*  cos(0*^  -0o)=  1 


(4.9) 
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In  Figure  4.2,  the  complex  0  plane  is  shown  with  the  original  path  of  integration, 
the  saddle  point,  the  steepest  descent  path,  and  the  branch  cuts  associated  with  each  wave 
type.  Note  that  there  is  no  fluid  wave  branch  cut  because  of  the  choice  of  transformation; 
that  is,  a  =  ko  cos  0  and  is  no  longer  discontinuous.  The  branch  cuts  are  placed  for  the 
marine  sand  case  (Vp  >  Vq,  Vs  <  Vq,  V2  <  Vq).  By  Cauchy's  theorem  the  path  of 
integration  in  the  complex  plane  may  be  deformed  arbitrarily  without  changing  the  value  of 
the  integral  provided  that  the  end  points  are  not  deformed  and  no  singularities  are  crossed 
from  one  path  to  the  other.  Taking  the  first  term  of  the  Taylor  series  (Brekovskikh  (1980), 
Eq.  27.9)  we  have 


Pr(9o)=^Po“' 


L^kgS 


il/2 


(4.10) 


Since  the  Taylor  series  is  expanded  in  orders,  l/k^s,  Eq  (4.10)  is  a  high-frequency,  far- 
field  approximation,  and  is  valid  for  k^s  »1.  The  region  of  validity  is  also  restricted  by 
the  approximation  of  a  slowly  varying  9l(0o).  so  that  Eq  (4.10)  will  not  be  valid  near 
pseudo-critical  angles  or  angles  of  pseudo-intromission.  Writing  the  asymptotic 
representation  of  the  Hankel  function  as 


(4.11) 


the  expression  in  Eq  (4.10)  for  the  reflected  pressure  can  be  considered  to  represent  an 
outgoing  cylindrical  wave  at  the  image  source,  modified  by  the  plane  wave  reflection 
coefficient,  to  first  order. 


Figure  4.2  Complex  9  plane  foi  the  reflected  field  showing  the  original 
contour;  the  saddle  point  ;the  steepest  descent  path  (SDP);  and 
the  Type  I,  Type  11,  and  shear  wave  branch  cuts  (BCp,  BC2, 
and  BCs  respectively). 
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4. 1.1. 3  Reflected  Field  (Head  Waves) 

When  the  angle  of  incidence  becomes  such  that  the  steepest  descent  path  crosses 
over  a  branch  point,  the  path  starts  and  ends  on  different  Reimann  sheets.  Cauchy's 
theorem  for  deforming  paths  is  then  violated  since  at  one  end  point  the  original  end  point  is 
on  the  upper  sheet  and  the  steepest  descent  path  is  on  the  lower  sheet.  In  order  to  remedy 
this  situation  the  path  of  steepest  descent  may  be  deformed  as  in  Figure  4.3  following 
Brekovskikh  (1980)  where  the  path  around  the  branch  cut  gives  the  head  wave  contribution 
and  the  rest  of  the  path  yields  the  reflected  field  approximated  by  Eq  (4. 10). 

1.  Type  I  Waves 

Considering  just  the  contribution  from  the  head  wave  field  arising  from  the  Type  I 
wave  branch  point  we  have  that 

-<«  0^ 

(4.12) 

where  iRjM.  is  the  reflection  coefficient  on  the  left  edge  of  the  cut  where  the  real  part  of  (k^p 
-  is  positive  and  the  9ip.  is  the  reflection  coefficient  on  the  right  side  where  the  real 

part  of  the  root  is  negative.  Along  the  branch  cut,  Im{(kp2  -  kx2)l/2}  =  0.  Combining  the 
limits  of  integration 


Php(x,z>0)  = 


-tPoto 


4jt 


,2  0 


Php(x,z>0)  =  -i^  r (Sp.- 

— joo 


(4.13) 


This  expression  can  be  evaluated  by  numerical  integration  and  the  simplest  w?  v  of 
numerically  treating  the  contour  is  to  transform  the  equation  back  to  the  kx  plane  such  that 
the  path  of  integration  becomes  a  straight  line; 


Figure  4.3  Deformation  of  steepest  descent  path  around  the  Type  I  wave 
branch  cut.  The  dotted  portion  of  the  path  is  on  the  lower 
Reimann  sheet. 


85 


P,p(x,z  ^  0)  =  ^  y  ( V  -  -  lcir'"e-<'‘.«dk. 

-too+k. 


(4.14) 


Thus  far,  the  generation  of  head  waves  for  the  case  where  the  Type  I  wave  velocity 
is  greater  than  the  upper  halfspace  velocity  has  been  considered  (Vp  >  Vq).  This  is  the 
wave  used  in  seismic  refraction  methods  to  estimate  sub-seafloor  layer  velocities  (see,  for 
example.  Clay  and  Medwin  (1977)).  We  now  consider  the  case  Vp  <  V©  which  results  in 
moving  ihe  branch  point  (and  cut)  from  the  fourth  to  the  second  quadrant  of  the  complex  0 
plane.  When  the  path  of  steepest  descent  crosses  the  branch  cut,  the  path  can  be  deformed 
as  before  around  the  cut,  and  the  head  wave  contribution  is  again  obtained  by  Eq  (4.13). 
For  this  case  the  head  wave  is  exponentially  damped  away  from  the  boundary.  The 
pseudo-critical  angle  is  complex  and  greater  than  Jt/2  when  Vp  <  Vq  and  is  obtained  by  Eq. 
(3.3). 

2.  Type  II  and  Shear  Waves 

There  is  also  a  Type  n  head  wave  for  both  cases:  V2  >  Vq,  and  V2  <  Vq.  The 
derivation  of  the  field  is  entirely  analogous  to  the  derivation  for  compressional  waves  so 
that 


Ph;(x.Z>0)= 


-9*2-) 


-i(ki-k’)"(z+h) 


(k2-k 


1) 


1/2 


-ik.x 


dk. 


(4.15) 


and  the  steepest  descent  path  first  crosses  the  Type  n  wave  branch  cut  at 
cosh  0lpj.2  cos  (0l^pj.2-0o)=l.  v/here 


®pe2 


=  sin 


Rc|k2l 

Re{k„} 
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The  contribution  firom  the  shear  wave  is  likewise  given  for  both  Vj  >  Vq  and  Vs  <  Vq  as 


Phs(x,z>0)  = 


-tpoto 

47t 


-9is_)e 


(4.16) 


and  the  steepest  descent  path  first  crosses  the  shear  wave  branch  cut  at 

cosh  cosCeRpcj  -  0o)  =  1  where  0pcs  is  given  in  Eq  (3.7).  The  shear  and  Type  U  head 

waves  decay  exponentially  away  from  the  boundary  when  the  velocity  is  less  than  the 
upper  halfspace. 

4.1. 1.4  Reflected  Field  (Surface  Waves) 

When  the  angle  of  incidence  becomes  such  that  the  steepest  descent  path  crosses 
over  a  pole  (or  poles)  of  the  integrand  Cauchy’s  theorem  requires  that  the  sum  of  the 
residues  of  the  pole(s)  (times  2jt  i^be  added  to  the  field.  The  physical  meaning  of  this 
contribution  is  the  phenomenon  of  surface  waves.  The  poles  arise  from  the  denominator  of 
the  reflection  coefficient  The  surface  wave  field  can  be  written  for  a  single  pole  as: 

Psurf(9o)  = 

or  in  terms  of  Cartesian  coordinates 
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where 


v-l/2 


res{9i(e)}=ies{9i(kj)(k5-k5)' 
and  the  spatial  dependence  of  the  field  is  contained  within  the  square  brackets.  For  a 


simple  pole 


res  91(0)j  = 


^,,here^(0)=^ 
(0)  A^.(0) 


(4.19) 


and  the  prime  denotes  a  derivative  with  respect  to  the  argument.  There  is  no  loss  due  to 
geometrical  spreading  expressed  in  Eq  (4.17,  4.18)  since  the  solution  is  confined  to  one 
horizontal  dimension  (surface  waves  spread  in  two  horizontal  dimensions).  The  intrinsic 
absorption  in  the  horizontal  dimension  is  given  by  the  term  exp(kl  x)  where  k^  is  the 
imaginary  part  of  the  surface  wave  wavenumber. 

The  location  of  the  surface  wave  poles  as  a  function  of  frequency  for  sand,  silty 
clay,  and  sandstone  is  shown  in  Figures  4.4  -  4.6,  respectively.  For  each  sediment  type 
the  pole  is  located  close  to  90°  on  the  real  axis  and  the  position  on  the  imaginary  axis  is 
strongly  a  function  of  sediment  type.  From  Figure  4.2  we  see  that  the  path  of  steepest 
descent  will  cross  the  silty  clay  surface  wave  pole  even  for  modest  angles  of  incidence.  A 
larger  angle  of  incidence  is  required  for  the  sand  surface  wave  pole  to  be  crossed  and  yet 
larger  for  the  sandstone  surface  wave  pole. 


4.1.2  Field  in  Lower  Halfspace 

For  the  field  in  the  lower  halfspace  we  must  first  choose  the  quantity  of  interest. 
The  most  common  measuring  device  used  in  the  seafloor  is  a  geophone  (that  measures 
displacement)  and  recently  in  the  underwater  acoustics  community  several  experimenters 
have  employed  a  hydrophone  (that  measures  the  pressure  field). 

rue  pressure  field  in  the  interstitial  fluid  is  given  by  Eq.  (2.25).  Using  the 
expression  *^or  u  and  w  in  Eq.  (2.65)  and  Eq.  (2.66) 
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Figure  4.5  Pseudo*Scholte  wave  pole  location  in  the  Complex  Angle  Plane 
for  Silty  Clay. 
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Figure  4.6  Pseudo-Scholte  wave  pole  location  in  the  Complex  Angle  Plane 
for  Sandstone. 
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.  oo  p  -i(ah+kjx) 

p(x,z<0)=^  J  [k2(MGp-C)Tpe'“'''  +  ki(MG2-  - - - dk, 


(4.20) 


The  total  displacement  field  D  is  given  by  D  =  u  +  U  =  2u-w/p  or: 


7  |ai{k,Tpc‘“l’"(2-Gp/p)+k,T2e'“"^(2-G2  /|3)+asV“>"(2-Gs  /P)  }+ 

a3{“p  +  Op  /  P)+a2?2e‘“""(-2  +  Gl  /  P)+kxTse“'"(2  -G,  /  p)  }  ]  x 

-i(cth+k,x) 

- - 

a 

(4.21) 

The  integral  expressions  for  the  transmitted  pressure  and  displacement  field  will  be 
evaluated  in  like  manner  as  the  reflected  field  by  considering  the  contribution  of  each  wave 
type  separately  including  the  geometrical  ray  path,  the  various  head  waves,  and  the  surface 
waves. 

4. 1,2.1.  Transmitted  Field  (Geometrical  Ray  Path) 

In  order  to  obtain  the  first-order  contribution  to  the  geometrical  ray  path,  we 
proceed  (as  with  the  reflected  field)  to  transform  the  integral  to  the  complex  9  plane.  In  the 
integral  representation  for  the  transmitted  field  consider  the  two  square  roots,  (ko^  - 
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and  (kp2  -  in  the  exponential.  The  latter  square  root  suggests  a  complex  angle 

transformation  kx  =  kp  sin  0  .  This  transformation  (by  proper  choice  of  coordinate 
system)  will  lead  to  an  integrand  commensurate  with  that  for  the  reflected  field  except  for  a 
square  root  in  the  denominator  (ko^  -  kp2  sin  0)^/2  giving  rise  to  a  pole  (collocated  with  a 
branch  point)  in  addition  to  the  poles  firom  the  transmission  coefficient  Tp(0).  In  order  to 
avoid  this  singularity,  the  transformation  kx  =  ko  sin  0  is  chosen.  The  next  step  is  the 
choice  of  a  suitable  coordinate  system,  because  the  coordinate  system  influences  the 
robustness  of  the  asymptotic  approximation.  Candel  and  Crance  (1981)  show  how  the 
choice  of  a  coordinate  system  centered  at  the  source  (as  was  done  for  the  reflected  field) 
yields  a  poor  approximation  to  the  transmitted  field  for  geometrical  ray  paths  near  and 
above  the  critical  angle.  They  suggest  rather  a  system  centered  at  a  virtual  source  depicted 
in  Figure  4.7  and  defined  by 


(4.22) 

x-Xp=SpSin0p 

Z-Zp=-SpCOS0p 


Sp- 


(x-Xpf+(z-Zpf 


1/2 


0 


p 


=  sm 


-I 


fw 

— sin0 

I  Vo 


o 


7 
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This  coordinate  system  will  be  employed  in  the  following.  Unlike  the  Candel  and  Crance 
example,  however,  the  lower  halfspace  of  interest  in  this  work  is  lossy  (i.e.  Vp  is  complex) 
such  tha^  angle  of  refraction,  9p,  and  the  reference  lengths  Sp,  Xp,  and  Zp  are  all  complex 
(see  Appendix  C). 

Rewriting  the  first  term  of  Eq.  (4.20)  now  with  the  transformation  kx  =  ko  sin  0 
and  the  coordinate  system  above,  we  have  that 


= 


47t 


71  . 


.-ik^hcose  -‘>toSine(spSinep+Xp) 
C  w 


d0 


(4.23) 


For  the  method  of  steepest  descent  let 

p=koSp 


f(0)=- 


— ^-sin^0 


cos0p+sin0  'in0p  > 


F(0) 


1/2 

=T  Zp^-iko(hcos0+XpSin0) 


(4.24) 
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The  saddle  point  is  then  defined  at 


f'(9)j0_0^  =O=cos0sin0p-sin0cos6 


— ^-sin^0 


-1/2 


(4.25) 


and  with  the  relation  (Snell's  law)  sin0p  =:^sin  0o  the  saddle  point  is 


0sp  —  Qq  •  (4.26) 

The  saddle  point  for  the  transmitted  field  in  this  coordinate  system  is  simply  the  angle  of 
incidence  corresponding  to  the  geometrical  ray  path.  There  is  a  subtlety,  however,  that 
bears  explication.  The  geometrical  ray  path  (as  shown  in  Figure  4.7)  connects  a  ray  in  the 
upper  medium  traveling  at  real  angle  0©  with  a  ray  in  the  lower  medium  at  real  angle  0ip. 
The  real  angle  of  transmission  for  a  lossy  medium  is  the  angle  associated  with  the  real  part 
of  the  wavevector  which  points  perpendicular  to  the  wavefronts  and  is  given  by 


f 

1/2 

01p=tan  * 

sin  00  / 

VoRe- 

< 

1 - ysin  ^00 

’o 

/^p 

• 

(4.27) 


while  the  angle  of  transmission,  0p,  associated  with  the  virtual  source  in  Eq.  (4.22)  is 
complex  and  given  by  Snell's  law. 
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Now  the  steepest  descent  path  is  given  by 
lm{f(0)}=lm{f(0sp)} 


(4.28) 


or 


f  \r2 


Re'll  ■^-sin^0 
<  P 

and  at  the  saddle  point. 


xl/2 


cos0p+sin0sin0p 


f(9o) 


=Re^-^^V„ 


(4.29) 


VqCOs  0p 


therefore 


(4.30) 


471*^^  ^  ^■ykoSpVp  COS0O 


Thus  from  Eq.  (4.20)  we  can  write 


(4.31) 


(4.32) 


p(x,z<0)=— 

47t 


where 
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:sin  ^ 


V 

— 2-sin0r 


y 


0o2=sin 


'' V 

-2- sin  09 

lV2 


'V 

y 


=kj(si-zjcos0i)+kQ(hcos0oi  +  Xjsin0oi)  ;  i  =  p,2 

Note  that  the  two  terms  in  Eq(4.32)  represent  two  ray  paths  at  two  distinct  angles.  Thus 
for  a  porous  medium,  the  first  order  predictions  consist  of  two  contributions  as  compared 
with  the  one  contribution  predicted  by  single  phase  theories. 

For  the  displacement  field* 


D— — ai{(2-Gp/p)kotan0op.|j^^^^  + 


(2-02  /P)k„un9„2,  ■;r^Pose2T2(e„2)e‘'''=  + 


(2-Gs/p)ksCOS0s 


COS0Q  —  /rt  1 _ 


koSsV^  cos  Bos 


33  (2-Gp/p)kpCos0 


2r.iV,  cos9p_  ,  ' 

P-\lkoSpVp  COS0OP  ' 


(2-G,/p)k2C0.92jgg:^?2(e„2)e-‘''- 

(2-G,/P)k„tan9<„  j|^cos9,2;(9o,)e-'’l''  }  ] 

^  koSsV 


(4.33) 
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4. 1.2.2  Transmitted  Field  (Head  Waves) 


As  in  the  reflected  field,  the  head  wave  contributions  come  from  an  integration 
around  the  branch  cuts  and  are  given  by 

(4.34) 


Mi 


MG„-C/_ 

Tn  e  P  -  Tp 


PoV2  I  Pp-'  -  >'P 

g-£(ah+k,x) 


MG2-C/-  _  -  \  ja2Z 

p„v| 


dk. 


a 


Ph2(x,z^0) 


-‘-^1 

-ioo+kj  ^ 


PoVp' 


g  -«(ah+k;ix) 


a 


dk,  (4.35) 


Phs(x.z<0)=  I 


47C 


-loo+k,  t- 


-Tp, + 

P„vJ  ^ 


MG2-C/-  - 

p„v| 


;  -'{ah+k,x) 


a 


dk. 


(4.36) 


The  sub-subscripts  refer  to  the  sign  of  the  vertical  wavenumbers  (ap,  ai,  as)  which  are 
real  on  their  respective  branch  cuts;  for  example  Tpp^  is  the  transmission  coefficient  for 

the  Type  I  wave  evaluated  on  the  left-hand  side  of  the  Type  I  wave  branch  cut 
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4. 1.2.3  Transmitted  Field  (Surface  Waves) 

As  in  the  reflected  field,  the  transmitted  field  due  to  a  surface  wave  can  be  written  in 
terms  of  the  residues  of  the  simple  poles  of  the  plane  wave  coefficients.  The  poles  for  the 
transmitted  field  are  identical  (co-located)  to  — jse  of  the  reflected  field.  For  the  pressure 
field  we  have: 


The  expressions  for  the  reflected  and  transmitted  field  are  in  the  form  of  Fourier 
integrals  such  that  the  Fast  Fourier  Transform  (FFT)  method  can  be  used  for  numerical 
solution  of  the  field  quantities.  For  example,  the  integral  in  Eq  (4.3)  can  be  approximated 
by  a  finite  Reimann  sum 


(4.38) 
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With  the  following  conditions  on  x,  kx  and  spatial  resolution  Ax: 

Xn,  =mAx  m  =  0,...,N-l 

kxn  =  nAkx  n  =  0,...,N-l 

Ax  =  27t/Ak,N  (4.39) 


we  have  that 


P(’‘ni.z)  = 


n=0 


4lt 


9l(nAkx) 


1/2 


(z+h) 


(ko-n^Ak^) 


-/27niin/N^ 


(4.40) 


This  expression  is  a  discrete  Fourier  transform  of  the  function  in  curly  brackets  and  may  be 
calculated  as  an  FFT. 

In  order  to  avoid  spectral  aliasing  the  sampling  frequency  ks  =  Akx  N  must  be  at 
least  twice  that  of  the  highest  frequency  to  be  transformed.  The  local  oscillation  frequency, 
however,  of  Eq  (4.3)  is 


-kx/^^(z  +  h)|= 


(z  +  h) 
T7I 


(4.41) 


which  tends  to  infinity  as  kx  approaches  k©,  so  that  some  aliasing  will  occur.  In  practice  ks 
is  chosen  a  few  times  greater  than  (k(z+h))  max. 

Another  problem  arises  due  to  the  denominator  in  Eq  (4.3)  which  goes  to  zero 
when  kx  =  ±ko-  This  integral  can  be  interpreted  as  Cauchy  principal  values  [Candel  and 
Oance  (1981)]  and  k©  should  divide  a  sampling  interval  in  equal  halves  such  that 
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ko  =  (J  +  1/2)  Akx 


(4.42) 


where  J  is  an  integer  less  than  N/4.  The  spatial  window  of  the  transformed  result  will  be 

L  =  NA  X  =  (J  +  1/2)  X  o  .  (4.43) 

In  the  following,  N  is  chosen  to  be  2048  and  J  is  200,  therefore  AX/X  ~  10  so  that  spatial 
wavenumbers  up  to  5  ko  will  not  be  aliased. 

4.3  Comparison  of  Two  Approaches 

In  this  section,  the  results  from  the  mixed  approach  are  compared  with  the  FFT 
results  in  both  the  upper  and  lower  halfspace. 


4.3.1.  Upper  Halfspace 

The  comparison  explicidy  stated  for  the  upper  halfspace  is 


3 


p  , —  —  ,  e  t(x(z+h) 

I  J 

j=l  — ioo+kj 

2n2res|2e-‘“<"'’>e-‘'‘-’‘ 


e"‘^**dkx  + 


j=l  —ioo+kj 


(4.44) 


where  indicates  inverse  Fourier  Transform,  Akx  is  the  sampling  interval  in  the  FFT, 


and  the  index  “j”  on  the  head  wave  summation  (second  term)  refers  to  the  three  wave  types: 
Type  I,  Type  II,  and  shear  wave.  The  plane  wave  reflection  coefficient  also  has  the 
subscript  “j”  in  this  term  to  indicate  that  along  the  path  of  integration  (i.e.  the  branch  cut) 
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the  vertical  wave  number  associated  with  each  wave  type  is  real,  and  is  positive  on  the  left 
side  of  the  cut  and  negative  on  the  right  side.  The  terms  on  the  RHS  of  Eq  (4.43) 
represent,  in  order,  the  geometrical  ray  path  contribution,  the  contributions  from  the  three 
head  waves  and  the  contributions  from  surface  wave(s).  The  quantity  (1+91)  was  chosen 
for  comparisons  because  it  is  functionally  less  complicated  than  91  and  was  consequently 
more  stable  in  the  FFT  calculations. 

Figure  4.8  shows  the  comparison  at  100  Hz  between  the  “exact”  (numerical 
integration)  solution,  the  FFT  result,  and  geometrical  ray  path  for  source  and  receiver 
height  of  50  meters  (about  3  wavelengths)  above  a  silty-clay  poro- viscoelastic  halfspace. 
For  this  geometry  and  sediment  type  the  contribution  of  surface  and  head  waves  are 
negligible.  In  addition,  the  first  order  geometrical  ray  approximation  quite  adequately 
predicts  the  field.  For  Figures  4.8  -  4.11,  the  horizontal  distance  is  expressed  in 
wavelengths  in  the  water. 

Figure  4.9  shows  a  comparison  at  100  Hz  of  the  FFT,  the  geometrical  ray  path,  and 
the  "exact"  solution  for  a  source  and  receiver  heights  of  20  meters  (about  1  wavelength) 
above  a  poro-viscoelastic  sand  interface.  For  this  geometry  and  sediment  type  the 
contribution  of  surface  and  head  waves  are  negligible.  The  geometrical  ray  approximation 
is  adequate  out  to  about  10  wavelengths.  At  this  range  the  angle  of  incidence  approaches 
the  critical  angle  and  the  assumption  that  91(0)  is  slowly  varying  breaks  down. 

Figure  4.10  shows  a  comparison  at  10  Hz  of  the  FFT,  the  surface  wave,  and  the 
two  head  waves  for  a  source  and  receiver  on  the  boundary  of  a  poro-viscoelastic  sand 
halfspace.  For  this  geometry  the  contribution  of  (1+91)  in  the  geometrical  ray  theory  limit 
is  zero.  The  difference  between  the  "exact"  result  and  the  summation  of  surface  and  head 
waves  is  due  to  the  fact  that  the  condition  of  a  slowly  varying  91(0)  is  violated  at  all  ranges 
since  the  geometrical  angle  of  incidence  is  90  degrees.  At  90  degrees  the  function  91(0)  is 
even  more  rapidly  changing  than  at  the  critical  angle.  Thus,  many  terms  in  the  expansion 
of  91(0)  are  required. 


Ab*tpr«««ur«|  (nrb) 
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Figure  4.8  -  Reflected  plus  image  source  field  at  100  Hz  calculated  by  FFT, 
numerical  integration,  and  asymptotic  expansion.  The  source  and  receiver 
are  50  m  above  a  poro-viscoelastic  silty  clay  halfspace. 


Abs|pr«ssure|  (arb) 


Figure  4.10  •  Reflected  plus  image  source  field  at  10  Hz  calculated  by  FFT 
and  numerical  integration.  The  source  and  receiver  are  on  the  boundary  of 
a  poro-viscoelastic  sand  halfspace. 
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4.3.2  Lower  Halfspace 

Figure  4.11  shows  a  comparison  at  100  Hz  of  the  FFT,  the  geometrical  ray  path, 
and  the  "exact"  solution.  Here  the  source  and  receiver  are  20  m  (about  1  wavelength) 
above  and  below  a  poro-viscoelastic  sand  interface,  respectively.  For  this  geometry  and 
sediment  type  the  contribution  of  surface  and  head  waves  are  negligible.  In  addition,  the 
geometrical  ray  approximation  quite  adequately  predicts  the  field. 

Predictions  for  the  transmitted  field  were  also  made  on  the  boundary  for  the  sand 
case  as  a  test  to  make  sure  that  the  boundary  conditions  were  met  numerically.  The 
pressure  field  predictions  for  the  transmitted  field  were  in  fact  virtually  identical  with  the 
sum  of  the  incident  and  reflected  field  results  on  the  interface. 

4.4  Summary  and  Conclusions 

In  this  chapter,  the  integral  expressions  for  the  field  from  a  line  source  above  a 
poro-viscoelastic  halfspace  were  evaluated;  by  a  decomposition  of  the  field  into  various 
wave  types  (the  "mixed"  approach),  and  by  FFT.  The  mixed  approach  results  were 
compared  with  the  exact  solution  and  the  FFT  results  in  order  to  determine  which  wave 
types  were  contributing  to  the  field  for  several  examples.  These  results  could  have  been 
compared  with  commensurate  quantities  using  the  single  phase  models,  however,  the 
differences  would  be  entirely  due  to  differences  in  the  plane  wave  reflection  and 
transmission  coefficients.  The  poro-viscoelastic  and  single-phase  models  of  plane  wave 
coefficients  were  previously  (in  Chapter  3)  compared  in  detail 


06  0« 
Ab«|pre(sure(  (arb) 
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Figure  4.11  -  Transmitted  fleld  at  100  Hz  calculated  by  FFT,  numerical 
integration,  and  asymptotic  expansion.  The  source  and  receiver  are 
respectively  20  m  above  and  below  a  poro*viscoelastic  sand  halfspace. 
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Chapter  5 

COMPARISONS  WITH  EXPERIMENTAL  DATA 

Recent  experiments  [Chotiros  (1989),  Satkowiak  (1988),  and  Muir  et  al.  (1979)]  of 
acoustic  interaction  with  the  seafloor  have  indicated  anomalous  behavior  of  high-frequency 
transmission  through  the  water- sediment  interface  when  the  surficial  sediments  are  granular 
(sands).  The  unexplained  results  include  anomalously  high  transmission  near  the  critical 
angle,  anomalously  low  phase  speeds  in  the  sediment  (-1000  m/s),  and  angles  of  refraction 
in  the  sediment  that  appear  not  to  obey  Snell’s  law.  The  anomalies  are  defined  as  lack  of 
agreement  with  classical,  single-phase  models  of  transmission  at  the  interface  between  two 
semi-infinite  halfspaces.  By  contrast,  the  experiments  indicated  good  agreement  with 
classical  single-phase  models  when  the  sediment  is  clay-like. 

One  of  the  anomalous  experimental  data  sets  will  be  compared  against  the 
predictions  for  transmission  through  a  poro-viscoelasdc  interface  to  determine  whether 
whether  the  Type  II  wave  plays  a  significant  role  in  the  transmission  process  and  is  the 
cause  of  the  anomalies.  The  acoustic  data  set  reported  by  Chotiros  (1989)  was  chosen  for 
analysis  in  this  chapter  because  the  accompanying,  comprehensive  sediment  geophysical 
and  geoacoustic  measurements  were  more  complete  than  for  any  of  the  other  data  sets. 

A  brief  description  of  the  experiment  [Chotiros  (1989)]  follows.  A  three- 
dimensional  array  of  12  hydrophones  was  buried  to  record  various  aspects  of  a  transmitted 
pulse.  The  hydrophones  were  arranged  in  two  orthogonal  ve^'ical  planes  on  two  sides  of  a 
cube  approximately  40  centimeters  in  length.  The  shallowest  probe  was  set  at  5  cm  below 
the  water-sediment  interface  at  the  intersection  of  the  two  planes.  The  source  was  mounted 
near  the  top  of  a  tower  (4.5  meters  above  the  seafloor)  that  could  be  moved  along  the 
seafloor  by  divers  to  vary  the  angle  of  incidence  of  the  impinging  sound  field  from  normal 
incidence  to  a  nominal  grazing  angle  of  22  degrees  (see  Figure  5.1).  The  frequencies  used 
in  the  experiment  were  5,  8, 15, 20,  30, 60,  and  80  kHz  and  the  associated  beamwidths  (at 
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the  3  dB  down  points)  varied  from  45  degrees  vertical  and  30  degrees  horizontal  at  the 
lowest  frequency  to  3  degrees  vertical  and  horizontal  at  the  highest  frequency.  The  raw 
acoustic  data  were  processed  for  signal  energy  densities  and  relative  arrival  times  using  a 
cross-correlation  technique  and  a  homomorphic  deconvolution  method;  plane  wave  front 
direction  and  group  velocity  using  a  least  squares  approach;  and  the  transmitted  sound 
pressure  level  (SPL)  which  required  essentially  no  signal  processing.  The  generated 
signals  included  various  pulse  shapes  and  lengths,  however,  for  the  data  of  interest  to  this 
chapter,  no  dependence  on  pulse  shape  or  length  was  found.  An  extensive  set  of 
geoacoustic  measurements  were  made  at  the  experimental  site  (the  particulars  of  the 
geoacoustic  data  are  discussed  in  detail  in  the  following  section,  section  5.1)  and  the  buried 
hydrophone  array  data  were  compared  with  the  SAFARI  [Schmidt  et  al.  (1985)]  model 
using  the  measured  sediment  geoacoustic  properties  as  inputs. 

The  receivers  buried  in  the  sediment  were  hydrophones  and  as  such  measured  the 
transmitted  pressure  field.  It  will  be  assumed  that  the  receiver  was  measuring  the  pressure 
field  in  the  interstitial  pore  fluid  (although  it  is  not  clear  that  the  fluctuations  of  the  solid 
sediment  frame  against  the  hydrophone  diaphragm  did  not  also  contribute  to  the  measured 
signal).  The  theoretical  transmitted  pressure  field  in  the  pore  fluid  was  derived  in  Chapter 
rv  and  is  given  by  Fxi  (4.20).  The  differences  between  the  formulation  in  this  work  (Eq 
(4.20))  and  the  experimental  data  just  described  include: 


Theoretical  Model 


EspsrimgntalPata 


Continuous  Wave 
Line  Source 


Pulse 

Beamed  source  with  beamwidth  a  function  of  frequency 


Ill 


Although  Eq  (4.20)  could  be  modified  to  take  into  account  the  effects  of  a  finite  pulse  and 
the  aperture,  the  processing  involved  to  make  predictions  directly  comparable  to  the 
experimental  data  seemed  fraught  with  difficulties,  the  resolution  of  which  seemed  beyond 
the  scope  of  this  work.  Thus,  the  predictions  of  Eq  (4.20)  will  be  used  to  examine  trends 
in  the  experimental  data  with  the  goal  of  determining  how  the  two-phase  nature  of  a 
granular  sediment  can  influence  transmissivity.  The  experimentally  measured  sound 
pressure  level  (SPL)  will  be  focused  on  since  this  was  a  direct  measurement.  The  other 
quantities  (wave  velocity,  transmission  coefficient,  refiraction  angle,  and  absorption)  were 
derived  quantities  (from  several  least-squares  and  deconvolution  techniques)  and  not  truly 
comparable  to  simple  notions  of  these  quantities  because  of  the  interference  of  the 
evanescent  field.  Comparison  of  these  derived  quantities  with  the  theory  herein  would 
require  using  the  ARL-UT  processing  programs  on  synthetically  generated  array  data.  This 
was  not  done. 

The  approach  taken  here  in  the  data  comparisons  is  to  first  derive  the  thirteen  Biot 
geophysical  parameters  for  the  sediment  at  the  experimental  site.  Some  of  the  parameters 
were  measured  directly,  some  of  the  parameters  can  be  obtained  by  empirical  or  semi- 
empirical  methods,  and  the  remaining  parameters  can  be  fitted  using  the  measured  acoustic 
data  (compressional  wave  speed  and  attenuation)  as  guides.  Biot  theory  is  then  used  to 
calculate  wave  speeds  and  attenuation  as  a  function  of  frequency  over  the  band  of  interest 
(the  experimental  data  were  obtained  at  frequencies  of  5  kHz  -  80  kHz).  These  results  will 
be  compared  with  the  geoacoustic  model  that  was  developed  based  on  the  core  data  and 
measurements  of  compressional  and  shear  wave  velocity  and  compressional  wave 
attenuation  across  several  cotes.  Finally,  the  anomalous  transmitted  SPL  data  near  the 
critical  angle  will  be  examined. 
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5.1  The  Sediment  Model 

In  this  section  the  13  Biot  input  parameters  are  obtained  for  the  sandy  sediment  at 
the  experimental  site. 

The  sediment  geoacoustic  properties  that  were  measured  included  the  comprcssional 
wave  velocity  and  attenuation  at  400  kHz  and  120  kHz,  and  shear  wave  velocity  at  around 
1  kHz.  The  measured  sediment  geophysical  properties  included  porosity,  permeability, 
and  grain  density  (K.  Briggs,  personal  cominunication).  The  grain  bulk  modulus  was 
taken  from  the  handbook  value  for  quartz  [Qark  (1966)].  The  properties  of  the  interstitial 
pore  fluid  were  detsimined  by  empirical  equations  based  on  the  temperature  and  salinity  at 
the  laboratory  conditions.  The  frame  shear  modulus  was  determined  from  the  product  of 
the  square  of  the  shear  velocity  and  the  density  which  is  nearly  an  exact  relation  since  the 
viscosity  of  seawater  is  small  (see  Eq(2.52)).  These  properties  (porosity,  permeability, 
grain  and  fluid  properties,  and  shear  modulus)  are  the  most  important  parameters  in  the 
Biot  theory  for  the  prediction  of  wave  speeds  and  attenuations;  the  remaining  parameters 
were  varied  within  reasonable  bounds  in  an  attempt  to  match  the  average  measured  wave 
compressional  velocity.  The  variability  in  predicted  compressional  velocity  with  variability 
in  these  parameters  was  quite  small,  with  predictions  varying  roughly  1825-1850  m/s  at 
400  kHz,  confirming  that  the  unknown  parameters  in  this  case  are  less  significant  than  the 
parameters  that  were  measured.  However,  no  combination  of  the  remaining  parameters 
would  produce  predictions  as  low  as  the  stated  1743  m/s  (sec  Table  5.1).  Thus,  the 
remaining  parameters  were  set  to  values  that  seemed  most  appropriate.  The  frame  bulk 
modulus  was  obtained  by  the  assumption  of  a  frame  Poisson’s  ratio  of  .2  [Hunter  (1961)]. 
The  pore  size  parameter  was  obtained  by  the  relation  [Hovem  et  al.  (1981)]: 


a  =  2(5  x/P)*'" 


(5.1) 
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The  logarithmic  decrements  for  both  the  firame  shear  and  bulk  moduli  were  chosen  to  be 
.05  which  is  roughly  in  the  middle  of  the  range  suggested  by  Stoll  (1980).  Finally,  the 
structure  factor  was  chosen  to  be  1.25  based  on  Stoll  (1974).  The  input  parameters  for 
Bic:  iheory  are  given  in  Table  5.2  and  the  predicted  velocities  and  attenuations  are  shown 
(along  with  the  measured  data)  in  Figures  5.2-5.7.  The  predicted  Type  I  wave  velocity  is 
somewhat  higher  than  that  measured,  the  predicted  Type  I  attenuation  is  substantially  lower 
than  that  measured,  and  the  predicted  shear  velocity  is  the  mean  of  the  measured  data  since 
the  frame  shear  modulus  was  chosen  from  the  mean  of  the  shear  velocity  data. 

It  is  believed  that  some,  if  not  all,  of  the  noted  differences  between  Biot  theory  and 
the  measured  data  are  attributable  to  scattering  in  the  core  measurements.  Theoretical  and 
experimental  evidence  [Fiona  and  Winkler  (1985)]  of  multiple  scattering  in  such  media 
indicates  that  compressional  velocity  decreases  monotonically  with  increasing  scatter,  and 
that  attenuation  increases  witii  increasing  scatter.  A  closer  look  at  the  measured  geoacoustic 
properties,  in  fact,  indicated  that  both  of  these  trends  exist  in  the  data.  As  evidence  of 
scattering  in  the  measurements,  consider  Figure  5.8  which  shows  measured  compressional 
velocity  at  120  kHz  and  400  kHz,  and  mean  grain  size  as  a  function  of  depth  in  the  core  at 
the  measurement  site.  First,  in  examiiung  just  the  frequency  dependence  of  the  two 
velocity  measurements,  the  negative  velocity  dispersion  is  clearly  evident;  that  is,  the  400 
kHz  velocity  measurement  is  consistently  lower  than  the  120  kHz  data.  Second,  scattering 
is  evident  from  the  middle  of  the  core.  Between  20  and  40  cm  there  was  a  layer  of  shell 
hash  which  is  evidenced  by  the  increase  in  mean  grain  size.  Note  that  the  measured 
velocity  at  400  kHz  decreases  sharply  with  increase  in  grain  size.^  For  granular 
sediments,  however,  compressional  velocity  is  theoretically  (and  empirically) 


^  Mean  grain  size  can  not  be  considered  to  be  perfectly  correlated  with  scattering  since  the  scatter  may  be 
caused  by  a  relatively  small  number  of  large  grains  that  may  not  be  well  indicated  in  a  mean  grain  size 
average. 


Table  5.1  ARL-UT  geoacoustic  model 


Property 

Units 

YalJi£ 

Fluid  density 

glcvr? 

1.02 

Fluid  speed 

m/s 

1536 

Sediment  density 

g/cm^ 

2.06 

Sediment  compressional 

m/s 

1743 

speed 

Sediment  compressional 

dB/m/kHz 

0.5 

attenuation 

Sediment  shear  speed 

m/s 

100 

Sediment  shear 

dB/na/kHz 

13.2 

attenuation 

Table  5.2  Geophysical  inputs  to  Biot  theory 


PropfirU 

Units 

Value 

Fluid  density,  (pg) 

g/cm^ 

1.024 

Fluid  bulk  modulus,  (Kf) 

dyne/cm^ 

2.396  lOlO 

Fluid  viscosity,  (ti) 

poise 

.0101 

Grain  density,  (Pj) 

g/cm^ 

2.66 

Grain  bulk  modulus,  (K^) 

dyne/cm^ 

3.6  10“ 

Porosity,  (p) 

.3673 

Permeability,  (k) 

cm2 

3.3  10-7 

Structure  factor,  (c) 

1.25 

Pore  size,  (a) 

cm 

4.2  10-3 

Frame  bulk  modulus,  (K^) 

dyne/cm2 

2.27  108+  i  3.6  106 

Frame  shear  bulk  modulus 
(li) 

dync/cm2 

1.7  108  + 1  2.7  106 

Type  I  Velocity  (m/i) 


Type  I  Attenuation  (dB/m/kHz) 


ISO  100  so 

Typo  II  Velocity  (m/e) 


200  190  100  90 

Typ*  n  Attenuation  (dB/m/kHi) 


Shear  Velocity  (m/e) 


Shear  Attenuation  (dB/m/kHs) 
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directly  proportional  to  grain  size.  This  negative  velocity  dispersion  in  the  data  strongly 
suggests  scattering.  Based  on  the  data  in  Figure  5.8,  it  is  clear  that  the  average  value  of 
corapressional  wave  velocity  of  1743  m/s  in  the  ARL-UT  model  (Table  5.1)  was 
contaminated  by  scatter.  Furthermore,  it  is  postulated  that  since  shell  fragments  were  noted 
throughout  the  sediment  column,  all  the  velocities  in  Figure  5.7  are  contaminated  by  scatter 
and  the  true  velocity  is  somewhat  higher,  perhaps  close  to  that  predicted  by  Biot  theory. 

In  Figure  5.9  the  attenuation  (in  dB/m/kHz)  is  plotted  against  depth  in  the  core  along 
with  the  mean  grain  size.  Consider  first  just  the  frequency  dependence  of  the  attenuation. 
The  consistent  increase  in  attenuation  with  frequency  is  again  strongly  indicative  that  the 
core  mea.<!’’rements  are  contaminated  by  scatter  (at  least  at  4(X)  kHz).  It  is  not  clear  without 
a  lower  frequency  measurement  whether  the  120  kHz  data  is  contaminated.  Second,  the 
large  (and  unrealistic)  measured  attenuations  toward  the  middle  of  the  core  are  almost 
certainly  highly  contaminated  by  scattering  losses. 

At  this  site,  group  velocities  in  ihe  granular  sediments  were  measured  at  about  1 1 00 
m/s.  It  was  suggested  by  Chotiros  (1989)  that  ihis  may  be  a  propagatory  wave  type  not 
predicted  by  classical  theory.  The  Biot  theory  in  fact  does  predict  another  wave  type  (in 
addition  to  the  classical  compressional  and  shear  velocity)  that  is  a  compressional  wave  and 
does  travel  at  lower  speeds  than  the  classical  compressional  wave  speed  (at  least  in  fluid- 
saturated  marine  sediments).  It  was  of  interest  therefore  to  calculate  its  velocity  to  examine 
the  possibility  of  whether  or  not  the  measured  wave  might  be  the  Type  H  Biot  wave.  At  80 
kHz  the  velocity  of  the  Type  n  was  predicted  to  be  roughly  100  m/s  and  no  reasonable 
combination  of  geophysical  parameters  could  force  predicted  velocity  as  high  as  that 
observed. 
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—  —  Mean  Grain  Size 
OOOOO  Velocity  400  kHz 
*^*^*-*^  Velocity  120  kHz 
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Figure  5.8  Measured  velocity  at  120  kHz  and  400  kHz  from  core  data 
taken  at  the  buried  array  site  with  mean  grain  size 
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—  —  Mean  Grain  Size 
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Figure  5.9  Measured  attenuation  at  120  kHz  and  400  kHz  from  core  data 
taken  at  the  buried  array  site  with  mean  grain  size 
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5.2  Sound  Pressure  Level  Comparisons 

The  experimental  transmitted  sound  pressure  levels  (SPL)  were  compared 
[Chotiros(l989)]  against  a  full  wave  model  (SAFARI)  that  considers  the  sediment  to  be  a 
solid.  The  geoacoustic  inputs  to  the  model  were  developed  from  the  core  data  and  are 
given  in  Table  5.1.  The  difference  between  the  measured  acoustic  transmissivity  and  the 
SAFARI  model  predictions  are  plotted  as  excess  SPL  in  Figure  5.10  for  two  tower 
positions  (TP);  TP  3  (corresponding  roughly  to  a  39  degree  grazing  angle)  and  TP  5 
(corresponding  roughly  to  a  22  degree  grazing  angle).  Positive  differences  in  Figure  5.10 
indicate  that  the  model  predictions  are  too  low.  Note  that  at  TP  3,  where  the  grazing  angle 
is  far  from  critical,  the  model  predictions  are  in  reasonable  agreement  with  the  data  except  at 
the  highest  frequency.  At  TP  5,  however,  the  measured  transmissivity  is  much  greater  than 
that  predicted  over  most  of  the  band.  The  standard  deviations  are  over  all  the  hydrophones 
in  the  array. 

It  is  the  purpose  of  this  section  to  determine  whether  the  Biot  theory  could  provide 
a  better  explanation  of  the  observed  transmissivity.  The  approach  taken  was  to  difference 
the  results  of  Eq(4.20)  using  the  geoacoustic  model  predicted  by  Biot  theory  from  Table 
5.2  (this  result  represents  the  Biot  theory  prediction)  with  the  results  of  the  first  term  of 
Eq(4.20)  using  the  ARL-UT  geoacoustic  model  (Table  5.1)  to  calculate  the  moduli,  the 
wavenumbers  and  the  transmission  coefficient  (which  represents  the  full-wave  solution 
where  the  sediment  is  considered  as  a  solid).  The  resulting  excess  SPL  then  can  be 
compared  with  Figure  5. 10. 

In  the  results  that  follow,  the  properties  of  the  water  are  chosen  to  be  representative 
of  in  situ  conditions  (density=  1.02  g/cm^  and  bulk  modulus=  2.4065  10^®  dyne/cm^ ). 
These  are  appropriate  for  comparisons  with  the  acoustic  data;  they  differ  somewhat  from 
the  fluid  properties  at  the  laboratory  conditions  in  Table  5.2,  and  are  consistent  with  the 
fluid  properties  in  the  ARL-UT  geoacoustic  model  (Table  5.1). 


SPL  ( 
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Comparisons  were  made  between  Eq  (4.20)  with  the  ARL-UT  geoacoustic  model 
and  Eq  (4.20)  using  Biot  theory.  For  the  comparisons,  a  receiver  depth  of  .4  m  was  used. 
The  ARL-UT  model-to-Biot  theory  comparison  results  are  shown  in  Figure  5.11  which 
compare  favorably  with  the  ARL-UT  model-to-data  comparison  in  Figure  5.10.  The 
excess  SPL  in  Figure  5.11  can  be  ascribed  to  one  or  a  combination  of  the  following  effects 
on  transmission  to  a  buried  receiver. 

excitation  of  Type  n  waves, 
surface  wave  contributions, 
head  wave  contributions, 

transmission  across  the  water-sediment  interface,  or 
attenuation  along  the  path  in  the  sediment 

In  order  to  isolate  which  effect(s)  caused  the  differences  observed  in  Figure  5.10, 
each  contribution  was  examined  separately.  The  second  term  of  Eq(4.20)  was  calculated 
to  determine  the  contribution  of  Type  II  waves  to  buried  receiver  and  was  found  to  be 
negligible  at  all  frequencies  at  both  tower  positions.  Next,  the  surface  and  head  wave 
contributions  were  calculated  from  Eq(4.37  and  4.34-4.36)  respectively;  the  pressure  field 
from  these  contributions  was  likewise  negligible. 

The  plane  wave  transmission  coefficient  was  calculated  for  Biot  theory  and  the 
ARL-UT  model.  The  magnitude  of  the  plane  wave  transmission  coefficients  are  shown  in 
Figure  5.12  where  at  a  given  angle  the  Biot  model  predicts  a  somewhat  lower 
transmissivity  due  to  conversion  of  energy  to  the  Type  n  wave.  The  conversion  is  stronger 
(hence  lower  transmissivity)  with  increasing  frequency.  Comparing  the  transmissivity 
between  the  Biot  and  the  ARL-UT  model  it  is  clear  that  the  differences  in  transmissivity 
between  the  two  models  accounts  for  only  a  few  dB  at  all  frequencies  and  at  both  sites. 
Thus,  the  transmission  differences  between  the  Biot  and  ARL-UT  model  are  not  enough  to 
account  for  the  excess  SPL. 


Crazinc  An|Ia  (deg) 


12  Plane  wave  transmission  coefficient  predicted  by  Biot  theory 
and  the  ARL*UT  geoacoustic  model 
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The  final  possibility  for  the  excess  SPL  in  Figure  5.1 1  is  differences  in  the  product 
of  sediment  path  length^  and  attenuation.  Indeed,  the  non-linear  attenuation  predicted  by 
Biot  theory  is  quite  different  from  the  ARL-UT  model  (see  Figure  5.13).  At  TP  3  the  path 
lengths  in  the  sediment  are  short  and  the  effect  of  the  attenuation  is  fairly  small.  Note, 
though,  that  the  effect  of  attenuation  increases  with  increasing  frequency  (since  attenuation, 
in  dB/m,  increases  with  increasing  frequency).  At  TP  5,  the  path  length  in  the  sediment  is 
longer  and  the  effect  of  the  different  attenuations  is  more  pronounced  (50  dB  excess  SPL  at 
80  kHz).  At  low  frequencies,  the  excess  SPL  in  Figure  5.1 1  shows  a  minimum  at  8  kHz 
similar  to  that  in  the  excess  data  SPL  (Figure  5.10).  The  cause  of  the  minimum  in  Figure 
5.11  is  the  crossing  of  the  Biot  non-linear  attenuation  with  the  linear  attenuation  of  the 
ARL-UT  model  between  8  kHz  and  15  kHz  as  shown  in  Figure  5.13.  That  is,  at  5  kHz 
and  8  kHz,  the  Biot  attenuation  is  greater  than  the  ARL-UT  attenuation  such  that  the  slope 
of  the  SPL  excess  is  negative;  at  15  kHz  and  above  the  Biot  attenuation  is  less  than  the 
ARL-UT  model  so  that  the  slope  of  the  SPL  excess  becomes  positive. 

5.3  Summary  and  Conclusions 

Data  were  examined  that  exhibited  anomalously  high  transmitted  SPL  (with  respect 
to  the  single-phase  model,  SAFARI)  in  order  to  determine  whether  the  anomalies  were  due 
to  the  Type  II  wave.  Calculations  indicated  that  the  Type  n  contributions  were  negligible 
and  were  not  the  cause  of  the  anomalies.  This  conclusion  rests  heavily  on  the  constructed 
geophysical  inputs  to  Biot  theory  for  the  sediment  in  question  (Table  5.2)  and  the  small  (~ 
100  m/s)  Type  II  wave  velocity  predicted  therefrom.  The  conclusion  is  robust,  though,  in 
that  many  of  the  geophysical  inputs  were  measured  and  that  varying  the  geophysical  inputs 


^  The  eflects  due  to  differences  in  the  sediment  path  lengths  predicted  by  Biot  theory  and  the  ARL-UT 
model  are  small  compared  with  the  effects  of  intrinsic  attenuation. 


Attanuation  (dB/m/kHs) 
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(within  physically  plausible  limits)  does  not  make  the  Type  n  wave  contributions  non- 
negligible.^ 

Much  of  the  anomalous  behavior  appeared  to  be  simply  the  result  of  an  incorrect 
compressional  wave  attenuation.  The  excess  SPL  at  TP  3  appears  to  be  entirely  explainable 
in  terms  of  the  compressional  wave  attenuation  (compare  Figures  5.10  and  5.11).  The 
excess  SPL  at  TP  5  appears  to  be  due  to  compressional  wave  attenuation  at  low  (5  kHz  and 
8  kHz)  and  high  (60  kHz  and  80  kHz)  frequencies.  The  two  mid-frequencies  (15  kHz  and 
30  kHz)  at  TP  5  are  not  well-explained  by  the  present  modeling  but  may  be  due  to  effects 
not  accounted  for  in  the  modeling  such  as  the  presence  of  the  shell  layer  (an  impedance 
boundary)  at  15  cm  (see  Figure  5.8).  Suffice  it  to  say,  however,  that  the  anomalies  are  not 
due  to  the  Type  II  wave. 

Thus,  in  large  measure,  the  failure  of  SAFARI  in  predicting  the  observed  levels  lay 
not  in  the  neglect  of  the  Type  n  waves  but  simply  in  the  incorrect  geoacoustic  model  which 
was  based  on  the  apparently  contaminated  core  data.  The  Biot  theory  provided  the  means 
to  check  the  core  data  based  upon  the  geophysical  inputs  that  arc,  in  general,  more  reliably 
measured. 


^  The  one  possible  exception  is  if  gases  were  present  in  the  interstitial  fluid.  Even  small  amounts  of  gas 
are  known  to  have  significant  effects  on  compressional  wave  velocity.  If  gas  were  present,  the  measured 
velocity  of  1100  m/s  may  be  velocity  of  the  Type  I  wave,  and  the  excess  SPL  in  Figure  S.IO  could  be 
explain^  by  simple  path  length  considerations. 
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Chapter  6 

SUMMARY  AND  CONCLUSIONS 

The  Biot  theory  of  propagation  in  a  porous  medium  provides  a  mathematical 
framework  for  studying  acoustic  interaction  with  the  seafloor.  The  theory  considers  the 
two-phase  porous  nature  of  marine  sediments  in  contrast  to  the  classical  models  of  wave 
propagation  in  the  seafloor  that  consider  marine  sediments  as  an  extended  single-phase 
fluid  or  solid.  A  boundary  value  problem  was  set  up  and  solved  for  a  line  source  in  a  fluid 
medium  above  a  poro-viscoelastic  halfspace.  Expressions  for  the  reflected  and  transmitted 
field  were  given  in  integral  form  and  asymptotic  expansions  in  the  high-frequency,  far-field 
limit.  A  set  of  simultaneous  equations  was  solved  to  give  plane  wave  reflection  and 
transmission  (Type  I ,  Type  n  and  shear  wave)  coefficients.  These  equations  also  yielded 
the  Scholte,  pseudo-Scholte,  and  pseudo-Rayleigh  wave  phase  velocities  and  attenuations. 

The  expressions  obtained  for  the  transmitted  field  in  the  pore  fluid  were  compared 
with  measured  data  from  an  array  of  hydrophones  buried  in  sand  and  insonified  at  high 
frequency.  This  data  set  had  shown  much  higher  transmissivities  near  the  critical  angle 
than  that  predicted  by  a  single-phase,  solid  model  (SAFARI).  The  Biot  theory  was  used 
both  to  develop  a  geoacoustic  model  (Type  I,  Type  n,  and  shear  wave  velocities  and 
attenuations,  and  density)  and  to  model  the  acoustic  transmissivity  of  the  sand.  The  Biot 
theory  showed  qualitative  agreement  with  the  measured  data.  Surprisingly,  this  was  not 
due  to  a  poro-viscoelastic  effect  (i.e.,  contributions  from  the  Type  n  wave)  as  speculated. 
In  fact  the  contributions  from  both  the  Type  II  wave  and  the  shear  wave  were  negligible. 
The  agreement  was  due  simply  to  the  realism  of  the  geoacoustic  model  (in  particular  the 
Type  I  wave  attenuation)  at  the  experimental  frequencies.  The  SAFARI  results  were  based 
on  measured  core  data  believed  to  be  contaminated  by  scatter. 

The  principal  objective  of  this  work  was  to  understand  how  the  porous  nature  of 
marine  sediments  influences  reflection  from  and  transmission  through  the  seafloor.  This 
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was  accomplished  by  comparing  plane  wave  coefficients  between  poro-viscoelastic  (using 
Biot  theory),  viscoelastic,  and  fluid  models  of  the  sediment.  The  differences  between  the 
geoacoustic  models  from  Biot  theory  and  from  the  single  phase  models  include:  the 
frequency  dependence  of  the  phase  velocities,  and  the  existence  of  the  compressional  wave 
of  the  second  kind.  The  effect  of  each  difference  on  the  plane  wave  coefficients  is 
summarized  in  the  following  paragraphs. 

The  frequency  dependence  of  the  velocities  and  attenuations  in  Biot  theory  were 
sufficient  to  cause  differences  in  the  plane  wave  coefficients  with  respect  to  those  predicted 
by  the  classical  single  phase  models.  For  the  sand  case  the  critical  frequency  was  low 
enough  to  be  in  the  frequency  range  of  interest.  Therefore,  the  velocity  and  attenuation  of 
each  wave  type  (especially  the  Type  I  wave  which  is  the  most  important  wave  influencing 
reflectivity  for  this  sediment  type)  varied  non-linearly  over  the  band  of  interest  leading  to 
the  frequency  dependent  plane  wave  coefficients.  For  the  silty  clay  case,  the  critical 
frequency  was  high  enough  that  the  velocity  of  the  Type  I  wave  (the  only  wave 
participating  strongly  in  the  reflection  and  transmission  process)  was  nearly  constant  with 
frequency  and  the  attenuation  was  slightly  non-linear  with  frequency.  The  reflection 
coefficient  was  independent  of  frequency  over  the  10  Hz  - 10  kHz  band.  For  the  sandstone 
case,  the  velocity  of  the  Type  I  and  shear  waves  were  only  weakly  dispersive.  From 
previous  results  this  suggests  that  the  plane  wave  coefficients  should  exhibit  only  a  weak 
frequency  dependence.  In  this  case,  however,  the  velocity  dispersion  of  the  Type  II 
significantly  influenced  the  plane  wave  coefficients  at  high  frequency. 

For  the  unconsolidated  sediments  considered  here  (silty  clay,  carbonate  sand,  and 
shelly  sand  of  Chapter  5)  the  excitation  of  the  Type  n  wave  was  negligible  at  all 
frequencies.  In  a  layered  sedimentary  sequence  (e.g.  turbidites  in  abyssal  plains)  with 
multiple  interfaces;  however,  the  excitation  of  the  Type  n  wave  may  provide  a  pon- 
negligible  loss  mechanism  not  contemplated  in  the  classical  single-phase  models. 
Moreover,  the  existence  of  such  layers  introduces  the  possibility  of  excited  resonances  of 
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the  Type  n  wave  in  the  layer.  In  any  case,  if  the  Type  n  wave  did  prove  to  be  important 
for  a  layered  medium,  it  would  only  be  important  at  high  frequencies  (with  respect  to  the 
critical  frequency)  where  the  Type  H  wave  is  propagatory.  The  losses  due  to  the  Type  II 
wave  were  clearly  evident  at  a  single  interface  for  the  consolidated  sediment  (sandstone)  at 
high  fr^uencies. 

Finally,  the  model-to-model  and  model-to-data  comparisons  provided  some 
unexpected  insight  into  the  propriety  of  the  common  use  of  core  data  (velocity  and 
attenuation)  in  developing  geoacoustic  models.  First,  core  data  measured  at  400  kHz 
cannot  in  the  general  case  be  extrapolated  in  frequency  by  assuming  a  linear  frequency 
dependence  on  attenuation.  Second,  the  core  data  for  sandy  sediments  may  be 
contaminated  by  scatter  (as  was  the  case  in  the  measurements  discussed  in  Chapter  V).  The 
Biot  theory  can  be  used  in  both  instances:  to  either  extrapolate  in  frequency  or  as  an  aid  in 
determining  the  reliability  of  the  core  data.  Examples  of  the  effect  of  extrapolation  from 
(synthetic)  core  data  at  400  kHz  were  given  for  plane  wave  reflectivity  and  surface  wave 
velocity  and  attenuation. 
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APPENDIX  A 

Definition  of  the  Complex  Elastic  Moduli 


The  complex  elastic  moduli  were  shown  by  Stoll  (1974)  to  depend  upon  the  bulk 
moduli  of  the  fluid,  solid,  and  frame  constituents  in  the  following  way: 


C  =  Ks  (Ks-Km)/(D-Kni) 


H  =  (Ks  -  Kni)2 /(D  -  Km)  +  Km+ 4  m/3 


M  =  Ks2  /(D-Km) 


D  =  Ks(l  +  P(Ks/Kf-l)  ) 

where  Ks  is  the  bulk  modulus  of  the  individual  grains,  Kf  is  the  bulk  modulus  of  the 
interstitial  fluid,  and  Km  is  the  complex  bulk  modulus  of  the  sediment  frame.  In  the 
general  case,  the  grain  bulk  modulus  can  be  complex  to  account  for  losses  within  the 
individual  grains;  however,  in  this  work  only  the  frame  bulk  modulus  and  the  frame  shear 
modulus,  M.  are  complex.  Losses  arising  from  a  non-ideal  interstitial  fluid  are  taken  into 
account  via  the  viscosity  rather  than  a  complex  modulus. 
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APPENDIX  B 
Numerical  Methods 

In  general,  the  computations  required  a  high  degree  of  precision.  This  is  primarily 
due  to  the  disparate  values  of  the  material  parameters  that  comprise  a  marine  sediment;  the 
difference  between  the  smallest  and  largest  parameter  typically  spans  15-20  orders  of 
magnitude  (e.g.,  see  Table  2.1).  This  resulted  in  numerical  instabilities,  for  example,  in 
the  solution  of  Eq  (2.70)  using  Cramer’s  rule.  Partial  pivotal  scaling  and  use  of  double 
precision  was  sufficient  render  the  matrix  calculations  stable. 

The  most  difficult  numerical  problem  encountered  was  the  solution  of  the  secular 
equation  which  required  finding  roots  of  a  transcendental  equation  in  the  complex  plane. 
Few  numerical  methods  for  solution  of  such  a  problem  in  fact  exist.  Muller's  (1956) 
method  was  attempted  but  was  non-convergent  even  at  double  precision  arithmetic.  Thus 
the  behavior  of  the  solution  surface  was  too  erratic  (under-resolved)  for  the  granularity  of 
the  number  system  of  the  dependent  variable.  Mathematica  ™  (Wolfram  (1988))  software, 
which  allows  arbitrary  precision,  was  then  employed  using  the  secant  method  extended  to 
the  complex  plane.  This  technique  was  successful  and  reasonably  efficient  (less  than  fifty 
iterations  were  usually  sufficient  for  convergence)  in  finding  the  roots  to  the  secular 
equation. 

The  numerical  implementation  of  the  theory  was  coded  in  Microsoft  Fortran  and 
executed  on  an  80386  personal  computer  with  an  80387  math  co-processor.  IMSL 
routines  were  employed  to  evaluate  the  Kelvin  function  and  its  first  derivative,  the  plane 
wave  coefficients,  and  the  FFT. 
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APPENDIX  C 

Derivation  of  the  location  of  the  virtual  source 

Candel  and  Qance  (1981)  derive  the  location  of  the  virtual  source  for  a  lossless 
two-fluid  boundary.  The  equations  for  a  porous  medium  can  be  rewritten  as: 


(C-1) 

where  the  subscript  i  =  p,  2  refer  to  the  Type  I  and  Type  II  wave,  respectively.  That  is,  for 
a  given  source/receiver  geometry,  the  angle  of  incidence  (0oi)  corresponding  to  the  steepest 
descent  geometrical  ray  path  will  be  different  for  each  wave  type.  Note  that  the  quantity 

’'o 


should  be  taken  such  that  the  imaginary  part  is  negative  (which  follows  from  the  radiation 
conditions  discussed  in  (Chapter  2).  Now  Eq.  (C-1)  together  with  Eq.  (4.22)  requires  that 
6oi  be  complex.  The  transcendental  equation  for  Ooi  can  be  written  as  : 
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z  sin0Qi  =  (h  tan0Qi  -  xjVj, 


Vi  .  2 

®oi 

*^0 


nl/2 


Xi 


(C-2) 


In  order  to  preserve  the  meaning  of  the  geometrical  ray  path,  the  following  changes  are 
made  to  the  virtual  source  position  as  derived  by  Candel  and  Crance  (1981).  The  angle  of 
the  ray  connecting  the  receiver  and  the  boundary  will  be  redefined  as  the  angle  of  the  planes 
of  constant  phase,  which  is  a  real  angle  and  given  by  Eq.  (4.27).  Eq.  (4.27)  can  be 
rewritten  to  depend  on  source/receiver  geometry  and  angle  of  incidence  as  follows : 


1/2 

z  sin0oi  =  (h  tan0oi  -  x)Vo  Re> 

’0 

X,- 

* 

(C-3) 


Note  that  the  angle  of  incidence  is  also  now  real.  For  the  results  in  Chapters  4  and  5,  Eq. 
(C-3)  is  first  solved  numerically  for  0oi,  the  x  coordinate  of  the  virtual  source  position  is 
then  given  by  Eq.  (C-1)  and  finally  the  z  coordinate  of  the  virtual  source  is  obtained  from 
Eq.  (4.22)  as 


(C-4) 


